감마선 이진성의 UV 라인 변동으로 보는 충돌풍 탐색

초록

본 논문은 회전 구동 펄서와 거대한 O/Be 별의 방출풍이 충돌해 형성되는 캐비티가 UV P Cygni 라인에 미치는 영향을 모델링한다. 캐비티의 개방각(η)과 관측선과의 각도(α), 그리고 궤도 거리(s)에 따라 라인 프로파일의 흡수·방출 구성이 어떻게 변하는지를 계산하고, LS 5039·PSR B1259‑63·LS I+61 303 세 시스템에 적용해 관측 가능한 변동성을 예측한다.

상세 분석

이 연구는 펄서 풍이 별풍을 밀어내어 형성되는 원뿔형 캐비티가 UV 공명선의 P Cygni 프로파일에 미치는 효과를 정량화한다. 기본 가정은 별풍이 β‑법칙(v(r)=v∞(1−R⋆/r)β)으로 가속되고, 펄서 풍은 거의 빛의 속도에 가까운 고에너지 입자 흐름으로 묘사된다. 두 풍의 램프 압력 균형점(Rs)은 p⋆=pₚ 식(1)으로 정의되며, 여기서 η=Ė/(Ṁ v(Rs) c) 가 펄서·별풍 압력비를 나타낸다. η<1이면 별풍이 우세해 펄서가 별 주변에 작은 캐비티를 만들고, η>1이면 캐비티가 별을 둘러싸는 넓은 원뿔이 된다. 캐비티의 반각 ψ는 수치 시뮬레이션 결과를 토대로 식(3)으로 근사한다.

라인 형성은 Lucy(1971) 방식을 확장해, 캐비티에 의해 차단된 영역을 제외하고 적분을 수행한다. 핵심 변수 fν(r)은 특정 속도면(Sν)에서 캐비티에 의해 가려진 면적 비율을 나타내며, 이는 η와 α(관측선과 캐비티 축 사이 각도) 및 Rs에 따라 달라진다. 프로파일 차이를 정량화하기 위해 A(η,Rs,α)=∫|P(v)−1|dv 로 정의한 면적을 사용하고, 단일 별 경우의 면적으로 정규화한 A₀를 평가한다.

파라미터 탐색 결과는 다음과 같다. (1) s가 작고 α가 작을수록, 즉 캐비티가 별에 가깝고 관측선과 정렬될수록 A₀가 크게 감소해 흡수 트로프가 얇아지고 방출 피크가 약해진다. (2) η가 1에 가까울 때 ψ≈90°가 되며, 캐비티가 별풍을 거의 절반 가려서 라인 형태가 가장 크게 변한다. (3) η가 매우 크면 캐비티가 별 전체를 둘러싸지만, 실제 물리적 제한(펄서 풍이 별 표면에 충돌) 때문에 η>~10⁴에서는 모델이 무의미해진다.

세 실제 시스템에 적용하면, LS 5039와 PSR B1259‑63는 현재 추정되는 η와 α가 각각 ψ≈90°·α≈90°(LS 5039)와 ψ≈30°·α≈180°(PSR B1259‑63)로, 캐비티가 관측선에 거의 가려지지 않아 A₀≈1에 가까워 변동이 거의 없을 것으로 예측한다. 이는 기존 HST 관측과 일치한다. 반면 LS I+61 303은 궤도 이심률과 경사각이 커서, superior conjunction(α≈180°)에서는 캐비티가 관측선을 크게 가리고, inferior conjunction(α≈0°)에서는 거의 가리지 않는다. 따라서 두 위상 사이에 A₀ 차이가 0.2–0.3 정도로 예상되어, UV P Cygni 라인의 흡수 깊이와 방출 강도가 눈에 띄게 변할 것이다.

관측 측면에서는 COS/HST의 분광 해상도(R≈3500)가 0.05–0.1 수준의 A₀ 변화를 감지할 수 있으므로, LS I+61 303에 대한 시계열 UV 관측은 펄서 풍-별풍 충돌 모델을 검증하는 강력한 시험대가 된다. 또한, 라인 이중성이나 이온화 구역(스트롬그레인) 효과를 무시한 단순 모델이므로, 실제 데이터와 비교할 때는 전리와 복사압, 비구형 캐비티 구조 등을 추가로 고려해야 한다.