LISA가 탐지할 은하계 조밀 이진 중력파 전경 예측

초록

본 연구는 최신 인구합성 모델을 이용해 은하계 내 백색왜성·중성자별·흑색별을 포함한 모든 조밀 이진의 중력파 신호를 예측한다. 시간·위치에 따라 변하는 금속성, 초신성 킥, 그리고 공통 외피(CE) 효율 α를 고려한 결과, LISA는 임무 시작 3개월 만에 전체 전경을 검출할 수 있으며, 4년 동안 2 000개 이상의 개별 이진을 식별할 수 있다. α 값에 따라 저질량(𝓜 < 1 M⊙)과 고질량(𝓜 ≥ 1 M⊙) 이진 비율이 변하므로, LISA 관측을 통해 CE 효율을 제약할 가능성을 제시한다.

상세 분석

이 논문은 은하계 조밀 이진의 전체 인구를 포괄적으로 모델링함으로써 LISA가 측정할 수 있는 중력파 전경을 최초로 정량화한다. 초기 10¹⁰개의 제로-연령 주계열(ZAMS) 이진을 Kroupa IMF와 균일 질량비(q) 분포, Sana et al. (2012)의 로그-주기·이심률 분포를 사용해 생성하고, 전체 별 질량 6.1 × 10⁹ M⊙에 해당하도록 설정하였다. 실제 은하 질량에 맞추기 위해 결과를 10배 스케일링하였다.

은하 구조는 bulge, thin disc, thick disc 세 구성요소로 나누어, 각각의 별 형성 역사(SFH), 금속성‑연령‑위치 관계, 그리고 내부 반지름·높이 스케일을 Wagg et al. (2022)의 반경-시간 의존 모델에 따라 할당하였다. 특히 thin disc의 반경 스케일이 시간에 따라 감소하도록 설정함으로써 내부-외부 성장(inside‑out growth)을 반영했다.

진화 단계에서는 최신 인구합성 코드 cosmic(v3.6.1)을 사용해 질량 손실, 질량이동, 조석, 초신성 킥, 중력파에 의한 궤도 감쇠 등을 모두 포함하였다. 초신성 킥은 Hobbs et al. (2005)의 Maxwellian σ=265 km s⁻¹를 기본으로, 전자‑캡처·초박형 초신성은 σ=20 km s⁻¹, 블랙홀은 프리백(fallback) 모델을 적용했다. 킥에 의해 발생한 중심질량 이동과 속도 변화를 galpy의 MilkyWayPotential2022로 3차원 궤도 적분하여 현재 위치를 재계산하였다. 이는 전경의 거리·위치 분포에 중요한 영향을 미치며, 특히 bulge 중심에서 킥된 시스템이 외부 디스크로 방출되는 현상을 재현한다.

CE 단계에서는 공통 외피 효율 α를 0.2, 0.5, 1.0 등 여러 값으로 변동시켜, 최종 이진의 궤도 주기·이심률·질량 분포에 미치는 영향을 조사했다. α가 클수록 에너지 전환 효율이 높아 외피가 더 쉽게 제거되고, 이로 인해 더 넓은 궤도와 높은 질량(특히 NS·BH 프라이머리) 이진이 살아남는다. 결과적으로 α가 증가하면 𝓜 < 1 M⊙인 저질량 이진 비율이 상승하고, 고질량 이진(𝓜 ≥ 1 M⊙) 비율은 상대적으로 감소한다. 이는 LISA가 측정하는 전경 스펙트럼의 형태와 개별 검출 가능한 소스 수에 직접적인 변화를 초래한다.

전파 신호 계산은 각 이진의 chirp mass 𝓜와 현재 궤도 주파수 f를 이용해 특성 진폭 h_c(f)= (2 G^{5/3}/c⁴ d) 𝓜^{5/3} (π f)^{2/3}을 구하고, LISA의 감도 곡선과 비교해 신호‑소음비(SNR)를 산출했다. 전체 전경은 파워 스펙트럼 형태로 합산되며, 3개월 이내에 LISA의 측정 노이즈보다 크게 나타난다. 개별 검출 가능한 이진은 SNR > 7 기준을 적용했으며, 4년 임무 동안 α=0.5 기준으로 약 2 200개가 검출될 것으로 예측된다.

핵심 인사이트는 다음과 같다. (1) 금속성·위치‑시간 의존 SFH와 초신성 킥을 포함한 전체 은하 모델링이 전경 강도를 10‑20 % 정도 상승시킨다. (2) CE 효율 α는 전경 스펙트럼의 저주파(대질량)와 고주파(저질량) 비율을 조절해, 관측된 스펙트럼 형태만으로도 α를 역추정할 수 있다. (3) 전경 검출이 빠르게 이루어지므로, LISA 초기 데이터에서도 전경 통계와 개별 소스 카운트를 활용해 은하계 이진 진화 모델을 검증할 수 있다.