Strongly interacting matter with criticality induced by modified excluded volume in core-collapse supernova simulations

기술 요약: 수정된 배제 부피(Modified Excluded Volume)에 의해 유도된 임계성을 포함한 강하게 상호작용하는 물질의 핵심 붕괴 초신성 시뮬레이션

문제 제기
핵붕괴 초신성(CCSN)의 폭발 메커니즘은 천체물리학에서 여전히 가장 중요한 미해결 과제 중 하나이며, 이는 고밀도 상태 방정식(EOS)에 크게 의존한다. 특정 시나리오는 강한 1차 상전이(hadronic matter에서 deconfined quark matter로의 전이)가 성공적인 폭발을 유발할 수 있다고 가정한다. 그러나 이전 연구들은 깁스 조건(Gibbs conditions)을 통해 결합된 별개의 강입자 물질 및 쿼크 물질 모델을 사용하는 "2상 접근법(two-phase approach)"에 의존해 왔다. 이 접근법은 밀도 도약(density jump), 발생 밀도, 그리고 상 사이의 압력 기울기에 대한 불확실성을 초래하며, 유한한 밀도와 온도에서 임계점을 갖지 못하는 경우가 많다. 또한, 최근의 중성자별 관측 결과(질량 $\approx 2 M_\odot$, 반경 $\approx 9\text{--}13$ km)와 중력파 이벤트(GW170817)에 따른 제약 조건은 핵물리학 및 천체물리학적 데이터와 일치하는 EOS 모델을 필요로 한다.

방법론
저자들은 밀도 의존적 메존-뉴클리온 결합(특히 DD2 모델)을 기반으로 하는 관계론적 평균장(RMF) 프레임워크의 새로운 다목적 EOS 클래스인 DD2-MEV를 개발하였다. 전통적인 2상 구성 대신, 저자들은 수정된 배제 부피(Modified Excluded Volume, MEV) 접근법을 사용한다.

• MEV 메커니즘: 이 모델은 배제 부피의 기하학적 해석을 대신하여 유효 자유도의 변화를 도입한다. 매질 의존적 배제 부피 함수 $\Phi_N(T, n_B)$는 뉴클리온의 퇴화 인자(degeneracy factors)를 조절한다. 이 함수는 반데르발스(Van der Waals) 거동을 모사하도록 설계되었으며, 과임계 및 미달 임계 영역과 고온에서의 임계 종단점(critical endpoint)을 특징으로 한다.
• 열역학적 일관성: 구현 과정에는 결합의 밀도 의존성으로부터 발생하는 재배열 항(rearrangement terms)이 포함되어 열역학적 일관성을 보장한다.
• 비교 모델: 연구는 두 가지 버전의 DD2-MEV EOS를 비교한다:
  1. DD2-MEV (Gibbs): 깁스 상 평형 조건(일정한 압력 및 화학 퍼텐셜)을 강제함으로써 반데르발스 불안정 영역을 제거한 버전.
  2. DD2-MEV (Van der Waals): 반데르발스 방정식의 특성인 음(-)의 압력 기울기와 불안정 영역을 유지하는 연속 모델.
• 시뮬레이션: 저자들은 AGILE-BOLTZTRAN 코드를 사용하여 구형 대칭에서의 일반 상대론적 뉴트리노 복사 유체역학 시뮬레이션을 수행한다. 6종류의 볼츠만 수송 스킴과 $40 M_\odot$ 전구체(s40a28)를 활용한다.
• 분석: 뉴트리노 신호에 대한 상세 분석과 GREAT(일반 상대론적 별진동학) 코드를 이용한 중력파(GW) 모드 분석을 수행하여 기본($f$) 모드와 중력($g$) 모드를 추적한다.

주요 기여 및 결과

1. 깁스 구성의 실패:
   표준 상전이 구성을 강제하는 DD2-MEV (Gibbs) EOS를 사용한 시뮬레이션은 초신성 폭발을 일으키는 데 실패했다. 이 모델에서의 강입자 상과 쿼크 상 사이의 밀도 도약은 이전의 RDF 하이브리드 모델보다 현저히 작다. 결과적으로 원생 중성자별(PNS)은 초음속 붕괴 대신 완만한 유체역학적 재구성만을 겪는다. 시뮬레이션은 충격파 회복 없이 블랙홀 형성으로 진행된다.

2. 연속적인 반데르발스 거동을 통한 폭발:
   대조적으로, 연속적인 DD2-MEV (Van der Waals) EOS를 사용한 시뮬레이션은 성공적인 폭발을 유도한다. 상전이는 중심 밀도와 온도의 급격한 상승을 유발하여 초음속 PNS 붕괴를 일으킨다. 그러나 이전 모델들이 상 경계에서 즉각적으로 강력한 제2 충격파를 형성했던 것과 달리, 이 모델은 지연된 충격파 형성을 보인다. 아음속 섭동파가 외부로 전파되어 몇 밀리초 후에 정지된 반동 충격파(standing bounce shock)에 도달한 후 강력한 충격파로 변한다.

   • 이젝타(Ejecta): 이 폭발은 이전 보고들보다 더 큰 질량($\sim 0.2 M_\odot$)을 초기에 방출하며, 진단 에너지(diagnostic energy)는 10 베테(10 Bethe, $10^{52}$ erg)를 초과하지만, 상당한 낙하(fallback)로 인해 최종 이젝타 질량은 $\sim 1 M_\odot$로 감소한다.

3. 뚜렷한 뉴트리노 신호:

   • Gibbs 케이스: 밀리초 단위의 버스트가 관찰되지 않는다. 뉴트리노 신호는 PNS의 구조적 재조정에 따른 단기적인 변동만을 보여주며, 이후 블랙홀 형성으로 이어지는 전형적인 강착(accretion) 거동을 따른다.
   • Van der Waals 케이스: 버스트 형태의 뉴트리노 신호가 방출되지만, 기존의 백 모델(bag models)이나 RDF EOS에서 보고된 1~2 ms 버스트에 비해 **상당히 넓은 폭(5–10 ms)**을 가진다. 이러한 확장은 지연된 제2 충격파의 형성과 전파에 기인한다. 이 버스트는 두 단계의 광도 상승을 특징으로 한다: PNS 수축에 의한 느린 초기 상승, 그리고 충격파가 뉴트리노스피어(neutrinospheres)를 통과할 때 발생하는 급격한 스파이크.

4. 중력파 모드 분석:
   연구는 PNS 진동($f$ 및 $2g_1$ 모드)에 대한 모드 분석을 수행한다.

   • 표준 피팅의 붕괴: 저자들은 모드 주파수와 PNS 평균 밀도($\bar{\rho}$) 및 컴팩트도($\kappa$)를 연결하는 표준 피팅 관계가 쿼크 물질의 자유도가 나타난 이후 붕괴됨을 발견했다.
   • 구별되는 거동: 실패한 (Gibbs) 케이스에서는 구조적 변화로 인해 $f$-모드의 기울기가 급격히 변하고 $2g_1$ 모드가 빠르게 감소한다. 성공한 (Van der Waals) 케이스에서는 두 모드 모두 상전이 및 충격파 형성 동안 큰 변동을 보이는데, 이는 이들의 흥분이 단순한 밀도/컴팩트도 스케일링과는 다른 물리적 현상에 의해 지배됨을 나타낸다.

의의 및 주장
본 논문은 수정된 배제 부피 메커니즘을 통해 1차 상전이를 모사하는 연속적인 EOS를 활용함으로써 전통적인 2상 접근법의 한계를 극복했다고 주장한다. 주요 의의는 상전이 구성의 세부 사항이 CCSN 폭발 메커니즘에 결정적임을 입증했다는 점에 있다.

• 연구는 "밀도 도약"이 핵심 파라미터임을 시사한다: 이 특정 MEV 모델의 깁스 구성처럼 작은 도약은 폭발을 방지하는 반면, 연속적인 반데르발스 거동은 폭발을 허용한다.
• 또한 뉴트리노 버스트 지속 시간과 중력파 모드 진화가 서로 다른 상전이 시나리오(실패 vs 성공, 연속적 vs 깁스 구성)를 구별할 수 있는 민감한 탐침임을 강조한다.
• 저자들은 이러한 결과가 기존의 QCD 기반 폭발 모델들과 큰 차이를 보이지만, 특정 파라미터화에 의존하고 있다는 점을 들어 신중하게 결론을 맺는다. 그들은 현재의 분석이 불안정 영역에서 음속을 영(0)으로 설정하여 가짜 해(spurious solutions)를 피했음을 언급하며, 다차원 시뮬레이션을 통해 이러한 결과들을 확인하고 직접적인 GW 신호를 계산할 필요가 있음을 강조한다.