Forming Double Neutron Stars using Detailed Binary Evolution Models with POSYDON: Comparison to the Galactic Systems

POSYDON 이진 인구 통계 합성 코드를 사용한 본 연구는 은하계 이중 중성자별이 특정 핵 정의 또는 높은 방출 효율 및 낮은 초신성 킥을 요구하는 두 가지 뚜렷한 공통 외피 하위 채널(Case B 및 Case C)을 통해 형성됨을 밝혀내며, 이를 통해 관측된 궤도 주기 분할을 설명하고 이 시스템들의 상세한 진화 물리학을 규명한다.

핵심

개요: 우주의 탐정 이야기

우주를 거대하고 혼란스러운 댄스 플로어라고 상상해 보세요. 대부분의 시간 동안 별들은 홀로 춤을 춥니다. 하지만 가끔 두 개의 거대한 별이 함께 태어나 서로 손을 잡기도 합니다. 수백만 년에 걸쳐 이들은 격렬하고 드라마틱한 춤을 추게 되며, 이는 두 가지 결말로 끝납니다. 영원히 멀어지거나, 아니면 서로 충돌하여 하나로 합쳐지는 것입니다.

이 논문은 **이중 중성자별(Double Neutron Stars, DNS)**에 관한 탐정 이야기입니다. 이들은 우주의 댄스 플로어에서 가장 무게감 있는 '헤비급' 선수들로, 매우 밀도가 높은 죽은 별(중성자별) 두 개가 서로 궤도를 돌고 있는 상태를 말합니다. 우리는 우리 은하인 은하수에서 이러한 쌍들을 약 25개 정도 알고 있습니다. 저자들은 POSYDON이라는 슈퍼컴퓨터 시뮬레이션을 사용하여, 이 쌍들이 정확히 어떻게 형성되는지, 그리고 왜 지금과 같은 모습을 띠게 되었는지를 밝혀냈습니다.

주요 발견: 같은 목적지로 향하는 두 가지 다른 경로

이 논문의 가장 놀라운 점은 이중 중성자별을 만드는 방법이 단 한 가지가 아니라는 것입니다. 저자들은 이 쌍들이 마치 같은 도시로 연결된 두 개의 서로 다른 고속도로처럼, **두 가지 뚜렷한 "하위 채널(sub-channels)"**을 통해 형성된다는 것을 발견했습니다.

그 차이는 드라마가 시작될 때 두 번째 별의 "나이"와 "성숙도"에 달려 있습니다.

1. "젊은" 경로 (헬륨 핵 채널):

   • 시나리오: 아직 성장할 시간이 충분하지 않은 젊은 별을 상상해 보세요. 이 별은 여전히 헬륨으로 이루어진 핵을 가지고 있습니다 (마치 아직 고등학교를 졸업하지 않은 십 대 소년처럼 말이죠).
   • 드라마: 이 젊은 별이 팽창하려고 할 때, 파트너(첫 번째 중성자별)에게 너무 가까워집니다. 그러면 이들은 **공통 외층(Common Envelope)**이라고 불리는 거대하고 지저식한 가스 구름 속에 뒤엉키게 됩니다.
   • 결과: 이 엉킴에서 살아남기 위해서는 가스 구름을 밀어낼 엄청난 에너지가 필요합니다. 만약 성공한다면, 이들은 매우 좁고 빠른 궤도에 놓이게 됩니다. 이들이 바로 결국 서로 충돌하여 중력파를 만들어낼(병합될) 쌍들입니다.

2. "늙은" 경로 (탄소-산소 핵 채널):

   • 시나리오: 이제 더 나이가 많고 성숙한 별을 상상해 보세요. 이 별은 연료를 다 태울 시간을 가졌으며, 탄소와 산소로 이루어진 핵을 가지고 있습니다 (마치 커리어를 마친 은퇴자처럼 말이죠).
   • 드라마: 이 늙은 별 역시 파트너와 함께 가스 구름 속에서 뒤엉키게 됩니다.
   • 결과: 늙은 별의 가스 구름은 더 느슨하고 밀어내기 쉽기 때문에, 시스템은 이 엉킴에서 훨씬 쉽게 살아남습니다. 하지만 이들은 훨씬 더 넓고 느린 궤도에 놓이게 됩니다. 이 쌍들은 우주의 수명 내에 서로 충돌하기에는 너무 멀리 떨어져 있습니다.

비유: 이것은 마치 붐비는 방 안에서 두 커플이 싸움을 말리려는 상황과 같습니다.

• **젊은 커플 (헬리움 핵)**은 좁고 답답한 공간에서 싸우고 있습니다. 군중을 밀어내는 데 많은 노력(에너지)이 필요하지만, 성공한다면 서로 아주 꽉 껴안은 상태가 됩니다.
• **늙은 커플 (탄소-산소 핵)**은 넓은 방에서 싸우고 있습니다. 군중을 밀어내기가 훨씬 쉽지만, 결과적으로 서로 멀리 떨어진 채 서 있게 됩니다.

"킥(Kick)" 문제: 왜 댄스 플로어가 미끄러운가?

별이 죽을 때, 별은 초신성 폭발을 일으킵니다. 이 폭발은 종종 새로운 중성자별에 마치 당구대의 수구(cue ball)를 치는 것과 같은 "킥(kick)"을 줍니다.

• 기존 이론: 과학자들은 예전에 이 킥이 (마치 당구 수구를 해머로 내리치는 것처럼) 매우 크고 무작위적이라고 생각했습니다.
• 새로운 발견: 저자들은 만약 킥이 그렇게 컸다면, 대부분의 쌍은 뿔뿔이 흩어져 이중 중성자별이 되지 못했을 것이라는 점을 발견했습니다.
• 해결책: 우리가 실제로 하늘에서 관측하는 모습과 일치하려면, "킥"은 훨씬 부드러워야 합니다 (마치 가볍게 툭 치는 것처럼 말이죠). 논문은 두 번째 별이 폭발하기 전에 몸무게의 대부분을 이미 잃었기 때문에, 두 번째 별이 받는 킥이 매우 부드러울 것이라고 제안합니다. 덕분에 이 쌍은 흩어지지 않고 결합을 유지할 수 있습니다.

첫 번째 별의 "재활용(Recycling)"

두 번째 별이 폭발하기 전, 첫 번째 중성자별(더 오래된 별)은 "메이크오버(외모 단장)"를 받습니다. 두 번째 별이 가스를 배출할 때, 첫 번째 별이 그 가스를 먹어치우는 것입니다. 이것은 마치 별을 더 빠르게 회전시키고 더 밝게 만드는 '우주적 다이어트 계획'과 같으며, 이를 통해서 "재활용된 펄서(recycled pulsar)"가 됩니다.

이 논문은 "젊은 경로(헬륨 핵)"가 첫 번째 별에게 더 좋은 식단을 제공하여 더 빠르게 회전시킨다는 것을 보여줍니다. 이는 우리가 하늘에서 보는 병합하는 쌍들(Sub-population i)이 매우 빠르게 회전하는 반면, 병합하지 않는 쌍들(Sub-population ii)은 더 느리게 회전하는 이유를 설명해 줍니다.

이상한 녀석들은 어떻게 된 걸까?

논문은 이 두 가지 경로 중 어디에도 속하지 않는 세 개의 이중 중성자별 그룹이 존재한다고 언급합니다. 이들은 설명하기 어려운 기이한 궤도를 가지고 있습니다. 저자들은 자신들의 모델이 이 세 개를 쉽게 설명할 수 없음을 인정하며, 이들이 (쌍으로 태어난 것이 아니라) 아마도 붐비는 성단 안에서 서로 충돌함으로써 형성되었을 가능성을 시사합니다.

결과 요요약

• 분리: 이중 중성자별은 두 가지 주요 경로를 통해 형성됩니다. 하나는 "젊은" 기증자(병합하는 쌍으로 이어짐)를 위한 경로이고, 다른 하나는 "늙은" 기증자(병합하지 않는 쌍으로 이어짐)를 위한 경로입니다.
• 에너지: "젊은" 경로가 작동하려면, 우주는 가스 구름을 던져버리는 데 매우 효율적이어야 합니다.
• 킥: 이 쌍들을 유지하기 위해서는 별을 만드는 폭발이 격렬하지 않고 부드러워야 합니다.
• 병합: 결국 충돌하여 중력파(GW170817 사건과 같은)를 만들어낼 쌍들은 거의 확실하게 "젊은 경로"에서 온 것입니다.

이 논문이 말하지 않는 것

• 다음 병합이 언제 일어날지 예측하지 않습니다.
• 언급된 세 개의 "이상한" 별을 포함하여 우주의 모든 기이한 별을 설명한다고 주장하지 않습니다.
• 별들이 시간이 지남에 따라 어떻게 "회전 속도가 줄어드는지(spin down)" 또는 어떻게 희미해지는지에 대한 모델을 포함하지 않으며, 이로 인해 현재 실세계의 관측값과 완벽하게 일치시키는 데 한계가 있습니다.

요약하자면, 이 논문은 상세한 컴퓨터 시뮬레이션을 사용하여 우주에는 이중 중성자별을 만드는 두 가지 뚜렷한 레시피가 있으며, 그중 오직 하나의 레시피만이 시공간의 구조를 흔들며 충돌하게 될 쌍들을 만들어낸다는 것을 보여줍니다.

기술 요약

기술 요약: POSYDON을 이용한 상세 이진 진화 모델 기반 이중 중성자별 형성 연구

문제 제기
이중 중성자별(DNS)은 거대 질량 이진 진화, 초신성 메커로니즘, 그리고 중력파원들을 이해하는 데 있어 핵심적인 실험실 역할을 한다. 은하계 내에서 20여 개 이상의 DNS가 확인되었으나, 그 형성 채널은 여전히 복잡하며 논쟁의 여지가 있다. 기존의 인구 통계 합성(population synthesis) 연구들은 주로 피팅 공식이나 단일 별 궤적 라이브러리를 사용하는 신속 이진 인구 합성(rapid BPS) 코드를 활용해 왔다. 이러한 접근 방식은 특히 공통 외포층(common envelope, CE) 단계와 핵심 붕괴 초신성(supernova, SN) 과정에서의 자기 일관적인 이진 상호작용 처리가 부족하다. 또한, 관측된 은하 DNS 인구는 궤도 주기와 이심률에서 뚜렷한 분기(bifurcation)를 보이는데, 표준 모델들은 미세 조정 없이는 이를 재현하는 데 어려움을 겪는다. 본 논문은 더욱 엄격하고 자기 일관적인 프레임워크를 통해 DNS 형성을 이끄는 상세한 진화 단계들을 규명해야 할 필요성을 다룬다.

방법론
저자들은 MESA 기반의 1차원 코드를 바탕으로 한 상세한 단일 및 이진 별 진화 모델을 통합한 최첨단 이진 인구 합성 코드인 POSYDON(POpulation SYnthesis with Detailed binary-evolution simulatiONs)을 사용한다. 신속 BPS 코드와 달리, POSYDON은 사전 계산된 별 진화 그리드와 온더플라이(on-the-fly) 계산을 사용하여 자기 브레이킹, 중력파 복사, CE 진화, 초신성 킥(kick)을 포함한 복잡한 물리학을 처리한다.

본 연구는 25개의 서로 다른 모델 변형(표 2에 요약됨)에 걸쳐 $2 \times 10^6$ 개의 이진 인구를 시뮬레이션한다. 조작된 주요 변수는 다음과 같다:

1. 공통 외포층(CE) 처리: CE 방출 효율($\alpha_{CE}$), 핵-외포층 경계의 정의(수소 질량 분율 1%, 10%, 30%), 그리고 포스트-CE 질량 손실 처리(안정적 질량 전이 vs 윈드 손실 vs 단일 단계).
2. 초신성(SN) 모델: 핵심 붕괴 결과(Sukhbold et al. 2016; Patton & Sukhbold 2020; Fryer et al. 2012)를 구분하여 철 핵 붕괴(CCSN)와 전자 포획 초신성(ECSN)을 구별하는 다양한 처방.
3. 출생 킥(Natal Kick) 처방: 표준 Hobbs-Maxwellian 분포($\sigma \approx 265$ km/s)와 "물리 기반" 모델들의 비교:
   • AsymEj: 방출 질량에 따라 스케일링되는 비대칭 방출 메커니즘(Janka 2017)에서 유도된 킥.
   • Linear: 방출 질량에 선형적으로 스케일링되는 킥(Richards et al. 2023).
   • Low Maxwellian: 감소된 분산 모델($\sigma \approx 20$ km/s).

모델들은 태양 금속함량을 가정하며 일정한 별 형성사를 가정한다. 결과로 도출된 합성 인구는 궤도 주기($P_{orb}$)와 이심률($e$) 분포, 그리고 병합율을 중심으로 관측된 25개의 은하 DNS 샘플(표 1)과 질적으로 비교된다.

주요 기여 및 결과

• 형성 채널의 분기: 본 연구는 첫 번째 생성된 중성자별과의 로슈 로브 넘침(Roche Lobe Overflow, RLO) 시작 시점에서 기증자 별의 진화 상태에 의해 결정되는 CE 형성 채널의 근본적인 분기를 식별한다:

  • He 핵 채널 (Case B RLO): 기증자가 헬륨 핵을 가진 상태에서 질량 전이를 시작한다. 이 시스템들은 일반적으로 더 짧은 궤도 주기를 가지며 더 높은 외포층 결합 에너지를 갖는다.
  • C/O 핵 채널 (Case C RLO): 기증자가 RLO 전에 탄소-산소 핵을 발달시킬 때까지 더 진화한다. 이 시스템들은 느슨하게 결합된 외포층을 가지며 더 넓은 궤도 주기를 결과로 낸다.
  • 의의: 이 분기는 은하 DNS 인구 사이의 짧은 주기 병합 시스템(Sub-population i)과 넓은 주기 비병합 시스템(Sub-population ii) 사이의 분리를 자연스럽게 설명한다.

• 공통 외포층 물리학에 대한 제약: (병합하는 DNS를 생성하는) He 핵 채널의 형성은 CE 효율에 매우 민알감적이다. 모델들은 이 채널이 다음 조건 중 하나를 충족할 때만 생존함을 밝혀냈다:

  • 관대한 핵 정의를 사용할 때 (H-분율 $\simeq 30\%$), 또는
  • 높은 CE 방출 효율을 채택할 때 ($\alpha_{CE} \gtrsim 1.2$).
    이러한 조건 없이는 He-핵 시스템이 CE 내에서 병합되어 사라진다. C/O 핵 채널은 더 낮은 $\alpha_{CE}$ 값에서도 생존할 만큼 더 견고하다.

• 출생 킥 제약: 표준 Hobbs-Maxwellian 킥 분포($\sigma \approx 265$ km/s)는 첫 번째 초신성 단계에서 대부분의 이진계를 파괴하여 관측된 것보다 수 자릿수 낮은 병합율을 초래하므로 배제된다. 관측된 은하 DNS 인구는 낮은 속도의 킥 처방($\lesssim 50$ km/s)에 의해 가장 잘 재현된다. 물리 기반의 AsymEj 및 Linear 모델은 박리된 전구체(stripped progenitors)에 대해 자연스럽게 낮은 킥을 생성하며, 궤도 분포와 국부 병합율 밀도($O(10) - O(100)$ Gpc$^{-3}$ yr$^{-1}$)를 모두 성공적으로 재현한다.

• 펄서 리사이클링 및 질량 강착: 연구에 따르면 펄서 리사이클링(첫 번째 생성된 NS의 회전 속도를 높이는 과정)은 CE 단계보다는 헬륨 별 동반자와의 안정적인 Case BB 질량 전이 단계 동안 주로 발생한다. He 핵 채널을 통해 형성된 시스템은 더 효율적인 리사이클링을 겪으며, 이는 관측된 짧은 회전 주기(Sub-population i)와 일치한다.

• 초신성 유형 및 전구체:

  • 첫 번째 초신성은 주로 철 핵 붕괴(CCSN)이다.
  • 두 번째 초신성은 종종 초박리(ultra-stripped) 전구체를 포함한다. 병합하는 DNS 전구체의 약 60%는 두 번째 SN 시점에 로슈 로브를 넘치고 있다.
  • 병합하는 DNS는 $\sim 0.2 M_\odot$ 의 중간 전 SN 외포층 질량을 가진 전구체로부터 태어난다. 이는 병합하는 DNS의 상당수가 Type Ib 또는 Type Ic 초신성으로 관측될 수 있음을 시사한다.

• 병합율 및 중력파원: 모델들은 He 핵 채널을 통해 형성된 DNS만이 허블 시간 내에 병합될 것이라고 예측한다. 따라서 GW170817 및 GW190425와 같은 중력파 사건은 이 채널에서 기원한 것으로 추론된다. 물리 기반 킥을 적용한 모델들의 예측된 국부 병합율은 LIGO-Virgo-Kaluera-KAGRA의 제약 및 은하 DNS 샘플과 일치한다.

의의 및 주장
본 논문은 DNS의 형성 역사를 해결하기 위해 상세한 진화 모델이 필수적이라고 주장한다. 피팅 공식을 넘어 자기 일관적인 그리드를 사용함으로써, 저자들은 관측된 은하 DNS의 다양성이 CE 단계 시작 시점의 기증자 별의 진화 상태로부터 오는 자연스러운 결과임을 입증한다.

저자들은 현재 POSYDON 내에 펄서 진화 모델(선택 효과 및 스핀다운을 고려하는 모델)이 부재하여 결론이 질적 비교에 기초하고 있음을 겸허히 언급한다. 또한 결과가 태양 금속함량에 특화되어 있음을 인정한다. 그럼에도 불구하고, 본 연구는 다음을 확립한다:

1. 형성 채널의 분기(He-core vs. C/O-core)는 관측된 궤도 주기 분리를 설명하는 이진 진화의 견고한 특징이다.
2. 효율적인 CE 방출(관대한 핵 정의 또는 높은 $\alpha_{CE}$를 통해)은 병합 인구를 형성하는 데 필요하다.
3. 낮은 출생 킥은 DNS 이진계의 생존을 위한 필수 조건이며, 이는 고속 킥 분포에서 벗어난 최근의 이론적 변화를 뒷받침한다.
4. 불안정한 Case BB 질량 전이는 이 상세 모델들에서 병합하는 DNS를 생성하지 않으며, 이는 초고속 병합 시스템 형성에 관한 기존 가설에 도전한다.

이 작업은 펄서 진화 모델이 POSYDON 프레임워크에 통합된 후 수행될 향후 정량적 베이지안 분석을 위한 토대를 제공한다.