NICER Perspective on TeV Blazar Mrk~421: X-ray Variability and Particle Acceleration

기술 요약: TeV 블레이저 Mrk 421에 대한 NICER의 관점: X선 가변성 및 입자 가속

문제 및 동기
블레이저, 특히 Markarian 421(Mrk 421)과 같은 고에너지 싱크로트론 피크(HSP) BL Lac 객체는 관측자의 시선 방향으로 정렬된 상대론적 제트에 의해 전자기 스펙트럼 전반에 걸쳐 상당한 가변성을 나타낸다. Mrk 421은 전형적인 TeV 방출체이자 가장 가까운 블레이저 중 하나($z=0.031$)이지만, X선 스펙트럼 진화와 급격한 가변성을 주도하는 물리적 메커니즘에 대한 상세한 이해는 여전히 활발한 연구 분야이다. 기존 연구들은 다양한 기기(예: Swift, XMM-Newton, NuSTAR)를 활용해 왔으나, 지속적인 기간 동안 NICER(Neutron Star Interior Composition Explorer)의 높은 감도와 높은 시간 분해능을 사용하여 Mrk 421의 X선 거동을 종합적으로 분석한 사례는 부족했다. 본 연구는 이 소스의 일간(intraday) 및 장기 가변성, 스펙트럼 곡률, 그리고 입자 가속 메커니즘을 규명하기 위한 스펙트럼 파라미터 간의 상관관계를 특징짓고자 하는 필요성을 다룬다.

방법론
저자들은 2022년부터 2024년 사이에 수집된 45개의 Mrk 421 NICER 관측 데이터를 분석하였다. 데이터 감소(data reduction)에는 HEASoft v6.34 및 XSPEC v12.14.1이 사용되었으며, 보정 최적화 및 배경 잡음 최소화를 위해 에너지 범위를 0.4–10 keV로 제한하였다.

• 타이밍 분석: 광도 곡선(light curves)은 60초 빈(bin) 크기로 추출되었다. 연구에서는 제곱평균제곱근(RMS) 가변성, 분율 가변성($F_{var}$)을 계산하였으며, 플럭스 분포를 조사하여 가우시안 또는 로그-노멀(log-normal) 프로파일을 따르는지 확인하였다.
• 스펙트럼 분석: 스펙트럼은 세 가지 모델(Power-Law (PL), Broken Power-Law (BPL), Log-Parabola (LP))을 사용하여 피팅되었다. 은하계 수소 칼럼 밀도는 $N_H = 0.0134 \times 10^{22} \text{ cm}^{-2}$로 고정되었다. 모델 선택은 감소된 카이제곱($\chi^2_r$)과 F-테스트를 사용하여 통계적 유의성을 결정하였다.
• 상관관계 연구: 저자들은 스펙트럼 파라미터(광자 지수 $\Gamma$ 또는 $\alpha$, 곡률 $\beta$, 싱크로트론 피크 에너지 $E_p$)와 플럭스 사이의 관계를 조사하였다. 스펙트럼 형태 변화를 평가하기 위해 경도비(Hardness Ratios, HR)를 계산하였다.
• 물리적 모델링: 스펙트럼 곡률을 해석하기 위해, 저자들은 관측된 $E_p$와 곡률 파라미터 사이의 반비례 관계를 재현하도록 로그-파라볼릭 전자 에너지 분포($N(\gamma) \propto \gamma^{-p-r \log_{10}(\gamma)}$)를 사용하여 싱크로트론 스펙트럼을 시뮬레이션하였다.

주요 결과

1. 가변성: Mrk 421은 2년의 기간 동안 강한 X선 가변성을 보였으며, 플럭스 수준은 약 28배(약 50에서 약 1380 counts s$^{-1}$까지) 변했다. 전체 관측 기간 동안의 총 분율 가변성은 약 48%였다. 일간 가변성 또한 유의미했으며, 몇 시간 단위의 시간 척도에서 변화가 발생했다.
2. 플럭스 분포: X선 계수율(count rates)의 분포는 다봉형(multimodal)이었으며, 저·중·고 플럭스 상태에 해당하는 세 개의 가우시안 성분의 중첩에 의해 가장 잘 피팅되었다. 이는 가변성이 단일 확률적 과정에 의해 구동되는 것이 아니라, 여러 개의 뚜렷한 방출 상태나 영역을 포함하고 있을 가능성이 높음을 시사한다.
3. 스펙트럼 모델링: 로그-파라볼라(LP) 모델이 45개 관측 중 42개에서 가장 정확한 설명을 제공하였으며, 이는 Power-Law 모델보다 성능이 현저히 우수했다. 나머지 3개의 관측은 Broken Power-Law에 의해 더 잘 설명되었다. LP 모델의 평균 스펙트럼 파라미터는 $\alpha \approx 2.32$ 및 $\beta \approx 0.32$였다.
4. 스펙트럼 경향:
   • Harder-when-brighter: 경도비(HR)와 계수율 사이에 강한 양의 상관관계($r=0.736$)가 발견되었으며, 이는 소스가 밝아질수록 스펙트럼이 단단해짐(hardens)을 나타낸다.
   • 파라미터 상관관계: 광자 지수 $\alpha$와 곡률 파라미터 $\beta$ 사이에 양의 상관관계($r=0.72$)가 관찰되었다. 반대로, $\beta$는 싱크로트론 피크 에너지 $E_p$와 음의 상관관계($r=-0.52$)를 보였고, $E_p$는 플럭스와 양의 상관관계($r=0.82$)를 보였다.
   • 피크 에너지 vs 곡률: 분석 결과, 싱크로트론 피크가 더 높은 에너지로 이동함에 따라(밝은 상태일 때), 스펙트럼 곡률($\beta$)이 감소함이 확인되었다.
5. 방출 영역 제약: 급격한 가변성 시간 척도($\sim40$분)와 싱크로트론 냉각을 가정했을 때, 자기장 세기는 $B \sim 0.5$ G로 추정되었다. 인과율 논거를 통해 방출 영역의 크기는 $R \lesssim 7.0 \times 10^{14}$ cm로 제한되었으며, 이는 X선이 중심 엔진 근처(수백 슈바르츠칠트 반지름 이내)의 조밀한 영역에서 기원함을 시사한다.

의의 및 주장
본 논문은 NICER 데이터를 기반으로 한 Mrk 421의 첫 번째 상세한 스펙트럼 및 타이밍 분석을 제시한다고 주장한다. 주요 의의는 HSP 블레이저의 X선 스펙트럼이 일관되게 로그-파라볼릭 모델로 가장 잘 설명된다는 점을 확인한 것이며, 이는 확률적 또는 에너지 의존적 입자 가속 시나리오를 뒷받침한다.

관측된 싱크로트론 피크 에너지($E_p$)와 곡률 파라미터($\beta$) 사이의 반비례 관계는 에너지 의존적 입자 가속(EDAP) 및 확률적 가속(SA)의 틀 안에서 해석된다. 이러한 모델들에서 입자 가속 확률은 입자 에너지가 증가함에 따라 감소하며, 이는 자연스럽게 로그-파라볼릭 전자 에너지 분포를 생성한다. 본 연구는 Mrk 421의 스펙트럼 진화와 가변성이 난류, 강한 자기장, 그리고 상대론적 유출을 포함하는 복잡한 비열적 과정에 의해 형성된다고 결론지으며, 여기서 스펙트럼의 곡률은 기저의 가속 메커니즘을 진단하는 도구 역할을 한다. 또한 다봉형 플럭스 분포는 제트가 단일 균질 구역보다는 여러 뚜렷한 방출 상태를 보유하고 있을 가능성을 암시한다.