DAO 새로운 공개 비상대론적 반사 모델

초록

DAO는 XSTAR와 Feautrier 방법을 결합해 고온 교정된 십자산 단면과 정확한 콤프턴 재분배 함수를 적용한 비상대론적 디스크 반사 스펙트럼 코드를 공개한다. 임의의 조명 스펙트럼을 지원하며, reflionx와 xillver와의 정밀 비교를 통해 높은 에너지 영역에서 기존 모델의 근사 한계를 보완한다.

상세 분석

본 논문은 X‑ray 바이너리와 AGN 디스크의 비상대론적 반사 스펙트럼을 계산하기 위한 새로운 공개 코드 DAO를 제시한다. 핵심은 XSTAR(원자 과정, 이온화 평형, 열 평형)와 Feautrier 방법을 이용한 2차 방사전달 방정식 해석을 결합한 점이다. 기존 reflionx와 xillver는 각각 Fokker‑Planck(코맥플레크)와 Gaussian 재분배 함수를 사용해 콤프턴 산란을 근사했으나, 고에너지(> keV) 혹은 고온(keV 수준) 전자 플라즈마에서는 이러한 근사가 크게 틀린다. DAO는 전자 분포를 완전한 상대론적 Maxwellian으로 가정하고, 에너지 의존적 σ_cs를 고온 교정하여 정확한 Klein‑Nishina 효과를 포함한다. 또한 Garcia et al. (2020)에서 제시한 정확한 재분배 함수 R(E_f,E_i)를 그대로 구현해, 한 번 이상 산란된 광자들의 에너지 이동을 양자역학적으로 처리한다.

DAO는 평면‑평행, 일정 밀도, 고정된 컬럼 밀도(10^24 cm⁻²) 가정 하에 azimuthal symmetry를 유지한다. 격자 설정은 τ = 10⁻⁴–10, 200 깊이 포인트, 5000 에너지 포인트(0.1 eV–1 MeV), 10 각도 포인트(코사인 0.05–0.95)이며, 주요 반복(N_main = 15)과 방사전달 반복(N_RTE = 200)으로 수렴을 보장한다. 입력 파라미터는 이온화 파라미터 ξ, 수소 밀도 n_H, 철 풍부도 A_Fe, 초기 온도 T_i, 조명 스펙트럼 형태(전력법, 전력법+컷오프, nthComp, 블랙바디) 등이다. 특히 하부(디스크 중간면)에서 블랙바디 방사선을 추가할 수 있어 XRB의 뜨거운 디스크를 모델링할 수 있다.

코드 구조는 4단계 반복: (1) XSTAR를 이용해 j(E), α(E) 계산, (2) Feautrier 방법으로 u(z,µ,E)와 h(z,µ,E) 구함, (3) 콤프턴 재분배 함수를 이용해 J_c(E) 업데이트, (4) 열 평형을 맞추기 위해 온도 프로파일을 재조정한다. 이 과정을 N_main 번 반복해 최종 방출 스펙트럼과 깊이별 이온 분포를 산출한다.

벤치마크 결과는 DAO가 reflionx와 xillver와 전반적으로 일치하지만, 특히 30 keV 이상에서 기존 모델이 과소평가하거나 과대평가하는 경향을 보인다. 정확한 재분배 함수를 적용했을 때, 고에너지 꼬리와 Fe Kα 라인의 폭·위치가 미세하게 변함을 확인했다. 또한, 조명 스펙트럼이 복잡할 경우(예: nthComp+블랙바디) DAO는 자유롭게 입력 파일을 받아 처리할 수 있어, 사용자 정의 조명 시나리오를 손쉽게 실험할 수 있다.

DAO는 완전 오픈소스(GitHub, Zenodo)이며, XSTAR와 동일한 원자 데이터베이스(ATDB 2024)를 사용한다. 따라서 향후 원자 데이터 업데이트가 자동으로 반영될 수 있다. 현재는 일정 밀도, 평면‑평행 가정에 머물지만, 코드 구조가 모듈화돼 있어 비등방성, 밀도 구배, 자기장 효과 등을 추가하기 용이하다.

요약하면, DAO는 기존 비상대론적 반사 모델의 핵심 한계인 콤프턴 산란 근사와 고온 교정 부재를 해소하고, 사용자 정의 조명 스펙트럼을 지원함으로써 X‑ray 반사 분석의 정확도와 유연성을 크게 향상시킨다.