Is the coexistence of strange quark stars and hadronic stars favored by astrophysical data? A Bayesian analysis

기술 요약: 이상 쿼크 별(Strange Quark Stars)과 강입자 별(Hadronic Stars)의 공존

문제 제기
미시적 구성 성분으로서의 컴팩트 별과 강하게 상호작용하는 밀집 물질의 상태 방정식(EOS)은 핵물리학 및 천체물리학의 핵심적인 미해결 과제로 남아 있다. 특정 긴장 관계가 존재하는데, 이는 중간 밀도에서의 "부드러운(soft)" EOS(일반 질량 천체의 작은 반경과 GW170817의 조석 변형성을 설명하기 위함)에 대한 관측 제약과, 고밀도에서의 "단단한(stiff)" EOS($2 M_\odot$를 초과하는 펄서들을 지탱하기 위함)에 대한 필요성 사이의 충돌이다. 더욱이, 강입자 물질에 하이퍼론(hyperons)과 $\Delta$ 공명(resonances)을 포함시키는 것은 일반적으로 EOS를 부드럽게 만들어, 최대 질량이 관측된 값보다 낮아지는 "하이퍼론 퍼즐(hyperon puzzle)"을 야기한다.

"두 가족(two-families, 2F)" 시나리오는 이를 해결하기 위한 방안을 제시한다: 두 개의 분리된 컴팩트 별 분과(branch)의 공존이다. 이 프레임워크에서 강입자 별(HSs)은 준안정 상태이며 작은 반경을 가진 분과를 점유하는 반면, 절대적으로 안정적인 이상 쿼크 물질(SQM)은 스스로 결합된 형태인 이상 쿼크 별(QSs)의 분과를 형성한다. 이 시나리오는 이론적으로 탐구되어 왔으나, 이전에 통합된 형태의 천체물리학적 및 실험실 데이터를 사용하여 표준적인 "한 가족(one-family, 1F)" 가설(모든 별이 강입자라는 가정)과 직접 비교하는 엄격한 베이지안 분석을 거치지는 못했다.

방법론
저자들은 1F와 2F 시나리오 모두에 대한 EOS 파라미터를 동시에 제약하기 위한 최초의 상세한 베이지안 분석을 수행한다.

• 모델:
  • 강입자 섹터: 핵자, 하이퍼론($\Lambda, \Sigma, \Xi$), 그리고 $\Delta$ 공명을 포함하는 비선형 상대론적 평균장(RMF) 모델을 사용한다. 이 모델은 포화 밀도에서의 핵 경험적 파라미터와 대칭 에너지 기울기($L_{sym}$)를 제어하기 위한 $\rho$ 메존의 밀도 의존적 결합을 활용한다.
  • 쿼크 섹터 (2F 전용): 절대적으로 안정적인 SQM을 기술하는 색-맛깔 결합(color-flavor-locked, CFL) 백 모델(bag-model) EOS를 사용한다. 이 섹터는 백 상수($B$), CFL 페어링 갭($\Delta_{CFL}$), 그리고 섭동적 QCD 보정($a_4$)을 포함한다.
• 데이터 및 가능도(Likelihoods): 분석에는 다음이 포함된다:
  • 질량-반경(M-R) 제약: NICER (PSR J0030+0451, PSR J0740+6620, PSR J0614–3329) 및 HESS J1731–347의 X선 관측으로부터 얻은 결합 사후 분포.
  • 질량 전용 제약: PSR J0952–0607 및 (잠재적 컴팩트 별로 취급된) GW190814의 이차 성분.
  • 중력파: GW170817의 전체 결합 사후 분포(질량 및 조석 변형성)를 사용하며, 2F 시나리오에서는 가족 레이블(HS-HS, HS-QS, QS-HS, QS-QS)에 대한 마진화(marginalization)를 적용한다.
  • 실험실/이론적 제약: 순수 중성자 물질에 대한 카이랄 유효장 이론($\chi$EFT) 밴드 및 대칭 핵 물질에 대한 중이온 충돌(HIC) 흐름 데이터.
• 통계적 프레임워크:
  • 본 연구는 사후 분포와 베이지안 증거($Z$)를 계산하기 위해 동적 중첩 샘플링(dynamic nested sampling, dynesty 패키지 사용)을 사용한다.
  • 모델 비교를 위해 베이즈 인자($B_{2F,1F} = Z_{2F}/Z_{1F}$)를 사용하며, 이때 서로 다른 파라미터 공간 차원(1F의 경우 11개, 2F의 경우 14개)을 고려한다.
  • 2F 시나리오에서 각 관측된 소스의 가족 할당은 잠재 변수(latent variable)로 취급되며, 구조적 안정성 및 천체물리학적 형성 경로에 기반한 사전 확률(priors)에 따라 마진화된다.

주요 기여

1. 최초의 베이지안 비교: 본 연구는 저밀도 핵 제약과 고질량 천체 관측을 모두 포함하는 포괄적인 데이터셋을 사용하여 1F 및 2F 시나리오를 엄격하게 비교하는 최초의 베이지안 분석을 제공한다.
2. 통합 프레임워크: 관측된 소스의 본질을 미리 가정하지 않고, 2F 가설 내에서 준안정 강입자 EOS와 절대적으로 안정적인 SQM EOS를 동시에 제약한다.
3. 긴장 해소: 본 분석은 2F 시나리오가 중간 밀도의 부드러움과 고밀도의 단단함 사이의 긴장을 두 개의 별 분과로 분리함으로써 자연스럽게 해결함을 보여준다.

결과

• 모델 선호도: 모든 데이터셋 조합에서 2F 시나리오가 1F 시나리오보다 선호된다. 전체 데이터셋(HIC 및 모든 천체물리학적 제약 포함)의 경우, 로그 베이즈 인자는 $\ln B_{2F,1F} = 11.9$로, 2F 시나리오에 대한 "결정적(decisive)" 지지를 나타낸다.
• EOS 특성:
  • 1F 시나리오: 작은 반경 제약과 $2 M_\odot$ 한계를 모두 만족시키기 위해, 1F EOS는 정교하게 균형 잡힌 타협을 요구한다. 이는 $\Delta$ 공점의 조기 출현($\sim 1.8 n_{sat}$)을 통해 달성되며, 이는 중간 밀도에서 EOS를 부드럽게 만들고 이후 고밀도에서 강력한 경화를 유도한다. 최대 질량은 $M_{max} \approx 2.05 M_\odot$로 제한된다.
  • 2F 시나리오: 긴장이 완화된다. 강입자 분과는 $\Delta$ 공명과 하이퍼론에 의해 지지되는 부드러운 상태를 유지하며, 최대 질량은 $M_{HS,max} \approx 1.54 M_\odot$, 반경은 $R_{1.4} \approx 11.1$ km이다. 무거운 펄서들은 더 단단한 QS 분과에 의해 지지되며, $M_{QS,max} \approx 2.33 M_\odot$에 도달한다.
• 파라미터 제약: 2F 사후 분포는 1F 사례에 비해 더 큰 디락 유효 질량($m^*/m$)과 덜 극단적인 $\Delta$ 결합 튜닝을 선호한다. 쿼크 섹터 파라미터($B, \Delta_{CFL}, a_4$)는 적절히 제약되며, 백 상수와 페어링 갭 사이에 양의 상관관계가 존재한다.
• GW170817 해석: 2F 시나리오에서 GW170817 이벤트는 HS-HS 병합일 가능성이 가장 높지만, 혼합 채널(HS-QS)도 상당한 기여(~35%의 가능도)를 하며, 최적의 EOS에 대해 QS-QS 채널은 기피된다.

의의 및 주장
본 논문은 현재의 천체물리학적 및 실험실 데이터가 단일한 순수 강입자 계열보다 두 개의 분리된 컴팩트 별 가족에 의해 더 자연스럽게 조직되어 있다고 주장한다. 2F 시나리오는 매우 컴팩트한 일반 질량 천체(HESS J1731–347 등)와 매우 무거운 펄서의 공존을 자연스럽게 설명하며, 하이퍼론 퍼즐을 효과적으로 제거한다.

저자들은 결론이 (하이퍼론 퍼텐셜, SQM 파라미터 범위, 소스 질량 가정 등의) 사전 확률 변화에 대해 견고함을 강조하지만, 제한된 데이터셋(HESS J1731–347, PSR J0614–3329 및 HIC 데이터 제외)만을 고려할 경우 통계적 선호도가 덜 결정적이라는 점을 언급한다.

본 연구는 현재의 비교에서 하이브리드 별(쿼크 핵을 가진 별)을 명시적으로 제외하였으며, 향-후 완전한 모델 선택을 위해서는 이들을 제3의 가설으로 포함하는 작업이 필요하다고 명시하였다. 저자들은 더 무거운 컴팩트 천체, 매우 컴팩트한 일반 질량 천체, 또는 동일한 반경을 가진 두 개의 뚜렷한 질량을 정밀하게 측정하는 것이 2F 시나리오를 지지하는 더 강력한 징후를 제공할 것이라고 결론짓는다.