Constraining the neutron star-black hole merger rate

우주가 거대하고 어두운 대양이라고 상상해 보세요. 중력파는 중성자별이나 블랙홀 같은 거대한 천체들이 서로 충돌할 때 발생하는 물결입니다. 수년 동안 과학자들은 LIGO, Virgo, KAGRA와 같은 거대한 귀라고 불리는 검출기를 사용하여 이 물결을 "듣기" 위해 노력해 왔습니다.

이 물결을 찾기 위해 그들은 "매칭 필터링(matched filtering)"이라는 방법을 사용합니다. 이것은 시끄러운 방 안에서 특정 노래를 찾는 것과 비슷하다고 생각하면 됩니다. 여러분에게는 이미 알고 있는 노래들의 재생 목록(이를 템플릿이라고 부릅니다)이 있고, 방 안의 소음과 여러분의 재생 목록을 비교하여 일치하는 것이 나타나는지 확인하는 방식입니다.

문제점: 재생 목록에 중요한 특징이 빠져 있었습니다
지금까지 과학자들의 재생 목록에는 커다란 사각지대가 있었습니다. 그들은 블랙홀과 중성자가 함께 춤을 출 때, 마치 피겨 스케이트 선수가 똑바로 서서 회전하는 것처럼 완벽하게 동기화되어 회전한다고 가정했습니다.

하지만 실제로는 이 우주의 무용수들이 종종 비틀거립니다. 만약 블랙홀이 궤도에 대해 이상한 각도로 회전하고 있다면, 전체 시스템은 **세차 운동(precession)**을 하게 됩니다(마치 쓰러지기 직전의 팽이처럼 비틀거리는 현상입니다). 기존의 재생 목록에는 이 "비틀거림"이 포함된 노래가 없었습니다. 따라서 만약 비틀거리는 쌍이 충돌한다면, 과학자들의 귀는 그 소리가 경직된 재생 목록과 일치하지 않기 때문에 그 신호를 완전히 놓칠 수도 있었습니다.

이 논문의 저자들은 블랙홀과 중성자가 매우 다른 질량을 가지고 있기 때문에, 이 비틀거림이 실제로 매우 흔하며 매우 뚜렷한 "소리"를 만들어낸다는 점을 깨달았습니다. 이를 무시함으로써, 우리는 발생하는 충돌의 최대 **85%**를 놓치고 있었을지도 모릅니다.

해결책: 더 똑똑해진 새로운 재생 목록
연구진은 이러한 "비틀거리는" 신호들을 처음으로 재생 목록에 포함하는 새로운 탐색 방법을 만들었습니다. 그들은 이 새로운 방법을 중력파 검출기의 세 번째 주요 관측 데이터에 적용하여 테스트했습니다.

그들이 발견한 내용은 다음과 같습니다:

초민감한 귀: 강하게 비틀거리는 시스템의 경우, 그들의 새로운 방법은 기존 방식보다 최대 100% 더 민감합니다. 이것은 마치 통조거 전화기에서 첨단 레이더로 업그레이드하는 것과 같습니다. 즉, 동일한 신호를 이전보다 두 배 더 멀리서도 들을 수 있습니다.
생각보다 적은 충돌 횟수: 이제 이 신호들을 훨씬 더 멀리서도 들을 수 있게 되었기 때문에, 그들은 우리가 듣고 있는 공간의 "부피"가 이전보다 훨씬 더 크다는 것을 깨달았습니다. 더 큰 부피의 공간을 관찰했음에도 여전히 적은 수의 충돌만 들린다면, 이는 실제로 우주에서 발생하는 충돌의 빈도가 이전에 계산했던 것보다 낮을 가능성이 높다는 것을 의미합니다. 구체적으로, 그들은 이러한 병합의 전체 비율이 이전 추정치보다 약 16% 더 작다는 것을 발견했습니다.

"비틀거리는" 하위 그룹
그들은 또한 "비틀거리는"(세차 운동을 하는) 쌍들을 구체적으로 살펴보았습니다. 초민감한 새로운 귀를 사용했음에도 불구하고, 그들은 데이터에서 확인된 비틀거리는 충돌을 단 하나도 발견하지 못했습니다. 이를 통해 그들은 엄격한 한계치를 설정할 수 있었습니다. 즉, 매년 10억 광년 입방 파섹당 발생하는 이러한 특정 비틀거리는 충돌은 79개를 넘지 않을 것입니다.

"거의 찾을 뻔한" 신호들
새로운 탐색은 실제 충돌로 확인되기에는 너무 미약한 네 개의 "경계선상의" 후보들을 포착했습니다. 흥고하게도, 이 네 가지 희미한 신호 모두 강한 비틀거림의 징후를 보였습니다. 그러나 과학자들은 신중합니다. 그들은 이것들이 실제 우주적 사건이라기보다는 지구에서 기원한 "정적"이나 노이즈일 가능성이 높다고 믿으며, 따라서 최종 수치에 포함시키지 않았습니다.

이것이 왜 중요한가
비틀거리는 회전을 포함하도록 "재생 목록"을 수정함으로써, 과학자들은 단순히 더 많은 신호를 찾는 것에 그치지 않고, 이러한 우주적 충돌이 얼마나 자주 발생하는지에 대한 더 정확한 그림을 얻고 있습니다. 이는 이러한 쌍들이 처음에 어떻게 형성되는지 이해하는 데 도움을 줍니다. 즉, 이들이 함께 평화롭게 진화한 별들로부터 태어난 것인지(보통은 거의 비틀거리지 않음), 아니면 밀집된 성단 속의 혼란스러운 별 무리에서 태어난 것인지(종종 비틀거림을 유발함)를 알 수 있게 해줍니다.

요약하자면, 그들은 더 나은 보청기를 만들었고, 우주가 생각보다 더 조용하다는 것을 깨달았으며, 비록 우리가 그 충돌 소리를 아직 듣지는 못했을지라도 우주의 무용수들이 예상보다 더 많이 비틀거리고 있다는 사실을 배웠습니다.

기술 요약: 중성자별-블랙홀 병합율 제약

문제 정의
LIGO–Virgo–KAGRA(LVK) 협업단이 수행하는 현재의 템플릿 기반 중력파(GW) 탐색은 사전 생성된 파형 뱅크에 대한 정합 필터링(matched filtering)에 의존한다. 이러한 기존 분석의 결정적인 한계는 일반 상대론적 스핀 유도 궤도 세차 운동(spin-induced orbital precession)을 무시한다는 점이다. 템플릿은 일반적으로 성분 스핀이 궤도 각운동량에 정렬되거나 반정렬된 시나리오로 제한된다.

이러한 제한은 질량 비가 커서 강한 스핀 세차 흔적을 보일 가능성이 높은 중성자별-블랙홀(NSBH) 쌍성계의 경우 문제가 된다. 만약 천체 물리학적 집단이 상당한 스핀 불일치(높은 유효 수직 스핀, $\chi_p$ )를 가진 쌍성계를 포함하고 있다면, 현재의 탐색 파이프라인은 이러한 사건의 상당 부분을 탐지하지 못할 수 있다. 결과적으로, 탐지된 NSBH 집단은 억제되거나 편향될 수 있으며, 이는 NSBH 병합율 및 근본적인 형성 경로(예: 고립 이진 진화 대 밀집 항성 환경에서의 역학적 포획)에 대한 부정확한 추론으로 이어질 수 있다.

방법론
저자들은 스핀 세차 운동을 명시적으로 포함하는 NSBH 쌍성계용 최초의 전용 템플릿 기반 정합 필터 탐색을 제시한다. 본 연구는 제3차 LIGO–Virgo–KAGRA 관측 주기(O3) 전체 데이터를 활용한다.

탐색 기법: 저자들은 불일치 스핀을 탐색에 포함하기 위해 새로운 방법(참조 문헌 [22]에서 제안됨)을 적용한다. 스핀 정렬을 가정하여 차원을 축소하는(두 개의 질량과 두 개의 스핀 크기로 구성된 4차원으로 문제를 축소하는) 표준 LVK 탐색과 달리, 이 "일반 스핀(generic-spin)" 탐색은 세차 항을 무시하지 않고 원형 궤도상의 NSBH 병합을 기술하는 데 필요한 전체 15차원 파라미터 공간을 고려한다.
민감도 추정: 탐색 민감도는 참조 집단 모델에서 추출된 시뮬레이션 신호를 실제 GW 검출기 노이즈에 주입하고 회복률을 계산함으로써 추정된다. 저자들은 근저에 깔린 NSBH 분포의 큰 불확실성으로 인해 불가지론적(agnostic) 집단 모델을 가정한다.
병합율 계산: 병합율 밀도는 추정된 탐색 민감도와 확신할 수 있는 GW 신호의 수를 결합하여 계산된다. 분석에는 천체 물리학적 비율에 대한 포아송 가능도(Poisson likelihood), 제프리스 사전 확률(Jeffreys prior)이 사용되며, 검출기 교정에서 발생하는 15%의 불확실성을 마진화(marginalize)한다.
비교: 저자들은 비교를 위한 베이스라인 역할을 할 수 있도록, 동일한 파라미터 공간을 다루면서 스핀을 정렬된 상태로 제한하는 축소된 버전의 표준 LVK 탐색을 생성했다. 본 연구의 모든 탐색은 가장 민감한 두 검출기인 LIGO-Hanford와 LIGO-Livingston에서 동시 검출된 신호로 제한되었다.

주요 결과

민감도 향상: 스핀 세차 운동을 포함하면 탐색 민감도가 크게 향상된다. 평균적으로 일반 스핀 탐색은 현재의 LVK 기술보다 44% 더 민감하다. 강한 세차 효과(유효 수직 스핀 $> 0.75$ )와 높은 총 질량( $> 14 M_\odot$ )을 보이는 계의 경우, 민감한 부피는 최대 100%까지 증가한다.
탐지 통계: 탐색을 통해 단 하나의 확신할 수 있는 탐지(FAR $< 10^{-2} \text{ yr}^{-1}$ )인 GW200115 042309를 식별했다. 베이지안 추론은 이 사건이 최소한의 세차 운동( $\chi_p \lesssim 0.4$ )을 가짐을 확인했다. Virgo 데이터를 제외했을 때 단일 검출기 후보였던 GW200105 162426는 회복되지 않았다.
병합율 제약:
- 세차 하위 집단: $\chi_p > 0.4$ 인 이론적 NSBH 쌍성계 하위 집단에 대해, 저자들은 $R_{90} = 79 \text{ Gpc}^{-3}\text{yr}^{-1}$ (90% 신뢰 구간)의 상한선을 설정한다. 이는 정렬된 스핀 파형만을 사용하는 탐색에 비해 상한선이 40% 감소한 수치이다.
- 전체 비율: 전체 파라미터 공간을 고려할 때, 국부 우주의 추정된 전체 NSBH 병합율은 $44^{+102}_{-37} \text{ Gpc}^{-3}\text{yr}^{-1}$ 이다. 이는 스핀 세차 운동을 무시하는 탐색 기술을 사용하여 얻은 비율( $52^{+130}_{-44} \text{ Gpc}^{-3}\text{yr}^{-1}$ )보다 평균적으로 16% 작다.
집단 분율: 세차하는 쌍성계의 비율(제로 관측을 가정할 때)을 예상되는 총 비율과 비교함으로써, 저자들은 우주의 NSBH 쌍성계 중 30% 이하만이 $\chi_p > 0.4$ 인 하위 집단으로서 스핀 세차 운동을 겪는다고 추정한다. 이 추정치는 정렬된 스핀 탐색으로부터 도출될 값보다 약 40% 낮다.
임계치 미만 후보: 탐색 과정에서 4개의 새로운 임계치 미만 후보(GW190421 203205, GW190916 134302, GW191225 091958, GW200310 030555)가 식별되었다. 이들은 모두 강한 스핀 세차 흔적을 보이지만, 천문학적 기원 확률( $p_{astro} < 1\%$ )이 낮아 지상 기원일 가능성이 높으며, 따라서 병합율 계산에서 제외되었다.

의의 및 주장
본 논문은 스핀 세차 운동을 고려함으로써 NSBH 탐색이 특히 등방성 스핀 분포를 가질 것으로 예상되는 역학적 경로를 통해 형성된 계에 대해 훨씬 더 민감해진다고 주장한다. 본 연구의 주요 의의는 국부 우주에서의 NSBH 병합율에 대해 더 엄격하고 정확한 제약을 가할 수 있다는 점이다.

저자들은 스핀 세차 운동을 포함함으로써 얻게 되는 증가된 민감한 부피를 고려하지 않았기 때문에 이전의 비율 추정치가 약 16% 과대평가되었을 가능성이 있다고 주장한다. 또한, 개선된 민감도는 NSBH 형성 메커니즘에 대한 더 엄격한 테스트를 가능하게 한다. 고도로 세차하는 후보가 관측되지 않았다는 사실을 통해, 저자들은 스핀 세차를 겪는 NSBH 집단의 분율에 대한 엄격한 상한선을 설정할 수 있다. 이러한 능력은 고립 이진 진화와 역학적 형성 경로를 구분하는 데 매우 중요하다. 본 논문은 이러한 결과가 지표적(indicative)이며 비율 추정의 불확실성을 모두 고려한 것은 아니지만, 스핀 세차를 무시하는 것이 집단 추론에 편향을 도입한다는 것을 보여준다고 명시한다.

유사한 논문