Binary neutron star populations across cosmic time: the impact of binary stellar evolution uncertainties

죽은 별들의 우주적 무도회

우주를 별들이 무용수처럼 춤을 추는 거대하고 혼란스러운 무도장이라고 상상해 보십시오. 대부분의 별은 홀로 춤을 추지만, 어떤 별들은 쌍으로 태어나 손을 맞잡고 공통의 중심을 돌며 회전합니다. 이 쌍들이 충분히 무거워지면, 결국 연료가 바닥나 붕괴하고 화려한 초신성 폭발을 일으키며 사라집니다. 만약 이들이 흩어지지 않고 폭발에서 살아남는다면, 도시 크기의 질량을 가진 중성자별이 됩니다. 이 별들은 너무 밀도가 높아서 티스푼 하나 분량의 질량이 10억 톤에 달할 정도입니다. 때때로, 이 두 중성자별은 서로를 꽉 붙잡은 채 궤도를 좁혀가며 점점 더 가까이 나선형으로 회전하다가 결국 서로 충돌합니다. 이 충돌은 우주에서 가장 격렬한 사건 중 하나로, 시공간에 '중력파'라는 파동을 보내며 금이나 백금 같은 무거운 원소들을 만들어냅니다.

오랫동안 우리는 우리 은하 내에서 이러한 "이중 중성자별" 커플을 단 몇 쌍밖에 알지 못했으며, 중력파 탐지 역사 전체를 통틀어 단 두 번의 충돌만을 목격했을 뿐입니다. 하지만 미래는 훨씬 더 시끄러워질 것입니다. 수천 건의 이러한 우주적 충돌을 매년 들을 수 있는 새로운 망원경과 탐지기들이 가동될 예정이기 때문입니다. 과학자들의 큰 질문은 이것입니다: 이 쌍들은 애초에 어떻게 형성되는가? 이는 별이 질량을 잃는 방식, 별들이 너무 가까워졌을 때 어떻게 상호작용하는지, 그리고 초신성 폭발로부터 발생하는 "킥(kick)"이 어떻게 그들을 흩어지게 만드는지를 포함하는 복잡한 퍼즐입니다. 이를 해결하기 위해 천문학자들은 우주의 시간 여행자 역할을 하는 슈퍼컴퓨터 시뮬레이션을 사용하여, 어떤 물리 법칙이 오늘날 우리가 보는 커플들을 만들어내는지 수백만 가지의 시나리오를 실행합니다.

위대한 우주 시뮬레이션

이 논문에서 저자 마티유 베네(Mathieu Venet)와 이리나 드보르킨(Irina Dvorkin)은 이 우주적 커플들을 대상으로 거대한 "만약에(what if)" 게임을 하기로 했습니다. 그들은 COSMIC이라는 강력한 컴퓨터 프로그램을 사용하여 우주의 역사 전반에 걸친 이중 성계의 삶을 시뮬레이션했습니다. 이것은 마치 물리 엔진을 미세하게 조정하여 이야기가 어떻게 변하는지 볼 수 있는 비디오 게임을 실행하는 것과 같습니다. 그들은 단순히 한 판의 게임을 돌린 것이 아니라, 별들이 질량을 교환하기로 결정하는 방식, 외층을 얼마나 효율적으로 벗겨낼 수 있는지, 새로 태어난 중성자별이 얼마나 강하게 걷차지는지, 그리고 그 킥이 정확히 언제 발생하는지라는 네 가지 특정 "조절 노브(knob)"를 변경하며 72가지의 서로 다른 버전을 실행했습니다.

그들은 환경이 결과에 어떤 영향을 미치는지 알아보기 위해 20개의 서로 다른 금속 함량 구간(별의 구성 성분 중 '무거운 원소'의 양을 나타내는 척도로, 매우 원시적인 금속 결핍 별부터 현대의 금속 풍부한 별까지 포함)에 걸쳐 시뮬레이션을 수행했습니다. 그들의 목표는 데이터 속에 이 별들의 진화 과정을 정확히 알려주는 보편적인 신호가 있는지, 아니면 불확실성이 너무 커서 한 시나리오를 다른 시나리오와 구별할 수 없는지를 밝혀내는 것이었습니다.

충돌로 가는 두 갈래 길

연구팀이 발견한 가장 놀라운 사실 중 중 하나는 이중 중성자별이 모두 동일한 레시피에서 나오지 않는다는 것입니다. 대신, 그들은 **이봉 분포(bimodal distribution)**를 발견했는데, 이는 별들이 중간에 "사막"을 두고 두 개의 뚜렷한 그룹으로 나뉜다는 것을 의미합니다.

작고 가벼운 재료를 사용하여 부드럽고 안정적인 과정을 거치는 케이크를 만드는 제과점을 상상해 보십시오. 다른 유형의 케이크는 거대하고 무거운 재료를 사용하여, 적절한 크기로 줄어들기 위해 혼란스럽고 무질서한 단계를 반드시 거쳐야 합니다.

• 저질량 그룹: 이 별들은 태양 질량의 약 6.5배에서 9배 정도로 가볍습니다. 이들은 대개 한 별이 다른 별을 삼켜버리는 혼란스러운 "공통 외층(common envelope)" 단계를 건너뛰며 부드러운 진화를 합니다. 이들은 보통 부드러운 "밀침(nudge)"을 주는 특정 유형의 초신성으로 폭발하여, 쌍이 흩어지지 않고 유지되도록 돕습니다.
• 고질량 그룹: 이들은 헤비급(태양 질량의 10~30배)입니다. 이들이 블랙홀이 아닌 중성자별이 되기 위해서는 반드시 혼란스러운 "공통 외층" 단계를 거쳐야 합니다. 이것은 마치 두 별이 서로 엉키는 우주적인 포옹과 같으며, 이 과정에서 거대 항성의 외층이 방출됩니다. 이 과정은 궤도를 극적으로 축소시키는데, 이는 매우 중요합니다. 왜냐하면 이 과정이 폭발 시 발생하는 격렬한 킥으로부터 쌍이 살아남을 수 있게 해주기 때문입니다.

그들 사이의 "사막"(태양 질량의 약 7~10배 구간)은 위험한 구역입니다. 이 범위의 별들은 저질량 그룹처럼 부드러운 밀침을 받기에는 너무 무겁고, 고질량 그룹처럼 혼란스러운 수축 과정을 견뎌내기에는 충분히 무겁지 않습니다. 이들은 보통 튕겨 나가 버려 이중 중성자별이 되지 못합니다.

불확실성의 조절 노브

저자들은 별의 유형은 일관되지만, 그 숫자는 물리 법칙의 노브를 어떻게 설정하느냐에 전적으로 달려 있다는 것을 발견했습니다. 가장 큰 두 가지 요인은 다음과 같습니다:

1. 킥 프리스크립션 (KICKFLAG): 이것은 새로 태어난 중성자별이 탄생할 때 얼마나 강하게 차여지는지를 제어합니다. 만약 킥이 너무 강하면(표준 설정), 대부분의 쌍은 흩어져 버립니다. 만약 킥이 더 약하거나 별의 질량에 따라 달라진다면(새로운 설정), 훨씬 더 많은 쌍이 살아남습니다. 이 단 하나의 변화가 이중 중성자별의 숫자를 **수 차례의 자릿수(orders of magnitude)**만큼 변화시킬 수 있습니다.
2. 공통 외층 효율 (ALPHA): 이것은 그 혼란스러운 "포옹" 단계 동안 별의 외층을 날려버리는 데 필요한 에너지가 얼마나 되는지를 제어합니다. 만약 효율이 높으면, 더 많은 쌍이 포옹에서 살아남아 함께 유지됩니다. 만약 낮으면, 그들은 합쳐져서 소멸합니다.

시뮬레이션 결과, 이러한 불확실성은 예측되는 이중 중성자별의 숫자를 100배 이상 변화시킬 수 있을 만큼 매우 컸습니다. 이는 현재 우리가 지금까지 목격한 몇 안 되는 병합 사례만으로는 어떤 버전의 물리학이 옳은지 확신할 수 없음을 의미합니다.

금속 함량의 반전

시뮬레이션의 혼란 속에서도, 저자들은 노브를 어떻게 설정하든 나타나는 몇 가지 "지문"을 발견했습니다. 이것들은 금속 함량이 변화함에 따라 형성율이 어떻게 변하는지에 대한 특정 패턴입니다.

• "건너뛰기(Skip)" 특징: 금속 함량이 매우 낮은 시기(초기 우주)에는 형성 효율이 급격히 떨어집니다. 이는 별들이 핵 연소 방식 때문에 특정 진화 단계(제1 거성 가지)를 갑자기 건너뛰기 시작하기 때문입니다. 이는 별이 지름길을 택함으로써 의도치 않게 이중 쌍 형성을 더 어렵게 만드는 것과 같습니다.
• "흔들림(Wobble)": 금속 함량 $10^{-3}$에서 $7.5 \times 10^{-3}$ 사이에서, 형성률은 단순히 올라가거나 내려가는 것이 아니라 비단조적(non-monotonic)으로 출렁입니다. 이는 별이 질량을 잃는 방식과 별들이 상호작용하는 방식 사이의 복잡한 상호작용의 직접적인 결과입니다.

결론

이 논문은 우리가 이 별들이 어떻게 형성될 수 있는지에 대한 좋은 지도는 가지고 있지만, 현재는 불확실성의 안개 속에 갇혀 있다고 결론짓습니다. 서로 다른 물리 법칙의 조합이 정확히 동일한 숫자의 병합을 만들어낼 수 있기 때문에, 어떤 규칙이 실제로 참인지 알기 어렵습니다. 이것은 마치 두 개의 서로 다른 열쇠가 같은 자물쇠를 여는 것과 같은 **퇴화(degeneracy)**의 문제입니다.

현재 중력파 탐지기는 이중 중성자별이 극도로 희귀한 최악의 시나리오들을 배제할 수 있는 수준에 있습니다. 이 결판을 내고 이 우주적 무도의 진정한 물리학을 밝혀내기 위해서는, 매년 최대 10만 건의 병합을 들을 수 있는 차세대 탐지기(2030년대 계획)를 기다려야 합니다. 이 방대한 데이터셋과 금속 함량이 시간에 따라 어떻게 변하는지에 대한 더 나은 측정이 결합될 때 비로소 우리는 노브를 미세하게 조정하여 이 무거운 커플들이 어떻게 모이게 되는지에 대한 진실된 이야기를 볼 수 있을 것입니다. 그때까지 우주는 시뮬레이션의 불확실성이라는 베일 뒤에 숨어 그 비밀을 간직할 것입니다.