[논문 리뷰] The 2.35 year itch of Cygnus OB2#9 III. X-ray and radio emission analysis based on three dimensional hydrodynamical modelling

E. R. Parkin, J. M. Pittard외 2명|arXiv (Cornell University)|2014. 06. 22.

High-pressure geophysics and materials참고 문헌 2인용 수 11

한 줄 요약

이 연구는 시리우스 OB2 #9 쌍성계의 3차원 적응형 메쉬 구조를 사용한 유체역학적 시뮬레이션을 통해 바람-바람 충돌 진동수에서 발생하는 X선 및 비열적 전파 방출을 모델링한다. X선 관측과 일치시키기 위해 낮은 전자-이온 온도 비율(~0.1)과 감소된 질량 손실률(~6.5–7.5×10⁻⁷ M☉ yr⁻¹)이 필요하다는 것을 발견하였으며, 전파 방출은 동기방출 방식이 주도적이며, 등가분포 가정이 성립하지 않는 에너지 밀도 비율을 요구한다. 이는 표준 모델을 초월한 복잡한 입자 가속 물리 현상을 시사한다.

ABSTRACT

X-ray and radio data recently acquired as part of a project to study Cyg OB2#9 are used to constrain physical models of the binary system, providing in-depth knowledge about the wind-wind collision and the thermal, and non-thermal, emission arising from the shocks. We use a three-dimensional, adaptive mesh refinement simulation (including wind acceleration, radiative cooling, and the orbital motion of the stars) to model the gas dynamics of the wind-wind collision. The simulation output is used as the basis for radiative transfer calculations considering the thermal X-ray emission and the thermal/non-thermal radio emission. To obtain good agreement with the X-ray observations, our initial mass-loss rate estimates require a down-shift by a factor of roughly 7.7 to $6.5 imes10^{-7}$ and $7.5 imes10^{-7}$ solar mass per year for the primary and secondary star, respectively. Furthermore, the low gas densities and high shock velocities in Cyg OB2#9 are suggestive of unequal electron and ion temperatures, and the X-ray analysis indicates that an (immediately post-shock) electron-ion temperature ratio of $\simeq 0.1$ is also required. The radio emission is dominated by (non-thermal) synchrotron emission. A parameter space exploration provides evidence against models assuming equipartition between magnetic and relativistic energy densities. However, fits of comparable quality can be attained with models having stark contrasts in the ratio of magnetic-to-relativistic energy densities. The radio models also reveal a subtle effect whereby inverse Compton cooling leads to an increase in emissivity as a result of the synchrotron characteristic frequency being significantly reduced. Finally, using the results of the radio analysis, we estimate the surface magnetic field strengths to be $\approx 0.3-52\;$G. (Abridged)

연구 동기 및 목표

• 다양한 파장에서의 관측을 통해 거대 쌍성 시리우스 OB2 #9의 바람-바람 충돌 영역 내 물리 조건을 제약하는 것.
• 충돌 가열 및 입자 가속의 맥락에서 X선 및 비열적 전파 방출의 기원을 조사하는 것.
• 등가분포 가정(예: 자기장 및 상대론적 에너지 밀도 간의 등가분포)이 이 시스템에서 성립하는지 테스트하는 것.
• 복사 냉각, 바람 가속 및 관측 각도가 관측된 광선 곡선 형성에 미치는 영향을 규명하는 것.
• 자기장 강도를 전파 방출 모델링을 통해 추정하고, 자석을 가진 별의 관측값과 비교하는 것.

제안 방법

• 두 별의 궤도 운동, 복사 냉각, 바람 가속을 포함한 3차원 적응형 메쉬 구조(AMR) 유체역학적 시뮬레이션을 사용.
• 유체역학적 시뮬레이션의 출력 결과를 바탕으로 복사 전달 계산을 수행하여 태양계 외부 X선 및 전파 방출 스펙트럼을 계산.
• 유체역학적 시뮬레이션에서 도출된 충돌 후 기체의 온도 및 밀도 분포를 기반으로 열적 X선 방출을 모델링.
• 상대론적 전자 에너지 분포 및 자기장 강도를 고려한 파rametric 주입 효율(ζ_rel, ζ_B)을 포함하여 동기방출을 통한 비열적 전파 방출을 시뮬레이션.
• 관측된 X선 및 전파 광선 곡선과 스펙트럼을 일치시키기 위해 ζ_rel, ζ_B 및 스펙트럼 지수 p의 파rameter 공간을 탐색.
• 전파 방출 모델링에서 자유-자유 흡수, 라진 효과, 역컴프턴 냉각 등의 감쇠 메커니즘을 고려.

실험 결과

연구 질문

• RQ1 시리우스 OB2 #9의 관측된 X선 스펙트럼을 재현하기 위해 충돌 후 영역에서 전자-이온 온도 비율은 어느 정도여야 하는가?
• RQ2 주성 및 보조성의 추정 질량 손실률은 이전 추정치와 어떻게 비교되며, X선 관측과 일치시키기 위해 어떤 조정이 필요한가?
• RQ3 관측된 비열적 전파 방출은 충돌 영역에서 자기장 및 상대론적 에너지 밀도 간의 등가분포와 일치하는가?
• RQ4 충돌 전 바람 속도 및 쌍성 간격의 변화가 궤도 단계 전반에 걸쳐 X선 및 전파 광선 곡선에 어떤 영향을 미치는가?
• RQ5 관측된 스펙트럼 지수 p ≪ 2를 초래할 수 있는 물리적 메커니즘은 무엇이며, 표준 확산 충돌 가속 예측과 비교해보면 어떻게 되는가?

주요 결과

• 관측된 X선 스펙트럼을 재현하기 위해 충돌 후 전자-이온 온도 비율이 약 0.1이 되어야 하며, 이는 전자 가열의 비효율성에 기인함을 시사한다.
• 이전 추정치 대비 질량 손실률이 크게 감소되어야 하며, 약 7.7배 감소된 값(주성: 6.5×10⁻⁷ M☉ yr⁻¹, 보조성: 7.5×10⁻⁷ M☉ yr⁻¹)이 필요하다.
• 전파 방출은 비열적 동기방출에 의해 주로 결정되며, 열적 방출은 무시할 수 없을 정도로 미미하여 이전의 전파 관측 결과를 확인한다.
• 자기장 및 상대론적 에너지 밀도 간의 등가분포를 가정한 모델은 스펙트럼 지수 p를 1에서 2 사이로 변화시켜도 관측 결과를 재현하지 못한다.
• 최적의 피팅 모델은 에너지 밀도 비율에서 극명한 대비를 보이며, ζ_rel ≫ ζ_B(예: ζ_rel = 0.15, ζ_B = 5×10⁻⁵, p = 1.6) 또는 ζ_rel ≪ ζ_B(예: ζ_rel = 3×10⁻⁴, ζ_B = 0.5, p = 1.2)의 형태를 취한다. 이는 표준 가속 물리 모델을 초월하는 복잡한 입자 가속 메커니즘을 시사한다.
• 표면 자기장 강도는 약 0.3–52 G로 추정되며, 강력한 자기장을 가진 단일 별의 일반적인 값보다 낮다. 관측된 p ≪ 2 값은 상류 입자가 충돌을 수정한다는 것과 일치하지 않으며, 다수의 약한 충돌에서의 재가속 메커니즘 등 대안적 가속 메커니즘을 시사한다.