Strange quark star I: the maximum gravitational mass and deformation of magnetized spinning model

기술 요약: 이상 쿼크 별 I: 자기화된 회전 모델의 최대 중력 질량 및 변형

문제 제기
본 논문은 극한의 조건(최대 $5 \times 10^{17}$ G의 강한 자기장 및 최대 1200 Hz의 빠른 회전) 하에 있는 컴팩트 천체에서 나타나는 이상 쿼크 별(Strange Quark Stars, SQS)의 구조적 파라미터를 다룬다. SQS(전체가 이상 쿼크 물질(SQM)로 구성된 가상의 컴팩트 천체)의 존재는 아직 확인되지 않았으나, 최근의 거대 컴팩트 천체 관측(예: PSR J0952-0607 및 GW190814 사건)은 표준 중성자별 모델에 도전 과제를 던지고 있다. 저자들은 밀도 의존적 상태 방정식(EOS)을 포함하는 SQS 모델이 이러한 관측 결과와 일치하는 중력 질량을 지지할 수 있는지 확인하고, 그에 따른 항성 변형 및 에너지론적 특성을 규명하고자 한다.

방법론
본 연구는 다각적인 이론적 및 수치적 접근 방식을 채택한다:

1. 상태 방정식 (EOS): 저자들은 가우시안 관계식에 의해 정의된 밀도 의존적 백 상수($B_{bag}(\rho)$)를 갖는 MIT 백 모델을 활용한다. 강한 내부 자기장을 고려하기 위해, 쿼크(업, 다운, 스트레인지)에 대한 페르미 관계식은 란다우 양자화(Landau quantization) 효과를 포함하도록 수정되었다. 모델은 제로 온도, 제로 전자 분율 근사를 가정하며, 전자의 기여는 $10^{-3}$ 차수이므로 무시한다.
2. 항성 구조 방정식: 3+1 형식론을 사용하여 축대칭 시공간에서 아인슈타인 장 방정식을 푼다. 에너지-운동량 텐서는 완전 유체, 자화, 그리고 순수 자기장의 기여를 포함한다. 이는 메트릭 함수들을 지배하는 4개의 타원형 편미분 방정식(PDE) 시스템을 생성한다.
3. 수치적 구현: 저자들은 스펙트럼 방법 기반의 C++ 코드인 LORENE 라이브러리를 사용하여 PDE를 해결한다. 계산 영역은 항성의 내부, 표면 근처의 얇은 껍질, 그리고 무한대로 뻗어나가는 외부의 세 영역으로 나뉜다.
4. 파라미터 공간: 회전 주파수는 0, 400, 800, 1200 Hz로, 중심 자기장 세기는 0에서 $5 \times 10^{17}$ G 범위로 설정하여 계산을 수행하였다. 각 구성에 대해 중심 엔탈피를 $0.01c^2$에서 $0.51c^2$까지 변화시키며 51개의 항성 모델을 산출하였다.
5. 검증: 수치 해의 정확도는 3차원 및 2차원 일반 상대론적 비리얼 항등식(GRV3 및 GRV2)을 사용하여 검증되었으며, 비자기/비회전 모델에서는 $\approx 10^{-5}$, 자기/고속 회전 모델에서는 $\approx 10^{-2}$의 오차를 보였다.

주요 기여 및 결과

• 최대 중력 질량 ($M_g^{max}$):

  • 모델은 비자기, 비회전 SQS에 대해 2.35 $M_\odot$의 최대 중력 질량을 예측한다.
  • 비회전 모델에서 자기장 세기를 높이면(최대 $5 \times 10^{17}$ G) $M_g^{max}$가 2.37 $M_\odot$로 약간 증가한다.
  • 회전은 최대 질량을 크게 향상시킨다. 가장 빠르게 회전하는 모델(1200 Hz)과 가장 강한 자기장($5 \times 10^{17}$ G)이 결합된 경우, 최대 중력 질량은 2.80 $M_\odot$에 도달한다.
  • 저자들은 자기 모멘트($\mu$)와 회전 주파수($f$)에 따른 $M_g^{max}$의 피팅 함수를 제공한다: $M_g^{max}(\mu, f)/M_g^{max}(0,0) \simeq (1 + a\mu^2)(1 + bf^c)$.

• 변형 파라미터 ($a$):

  • 변형 파라미터는 적도 반지름과 극 반지름의 비율로 정의된다 ($a = R_{eq}/R_{pol}$).
  • 연구 결과, 회전이 변형의 지배적인 요인이며, 자기장은 편평도(oblateness)에 미미한 영향을 미치는 것으로 나타났다.
  • 발견된 최대 변형 파라미터는 1.55이며, 이는 $f = 1200$ Hz 및 $B_c = 5 \times 10^{17}$ G인 모델에서 발생한다.
  • 변형에 대한 피팅 함수 또한 제공된다: $a(\mu, f)/a(0,0) \simeq (1 + \tilde{a}\mu^2)(1 + \tilde{b}f^{\tilde{c}})$.

• 에너지론 및 안정성:

  • 결합 에너지: 바리온당 총 결합 에너지($|E_{BE}|/A$)는 171 MeV에서 184 MeV 사이이며, 이는 SQS 모델의 안정성을 확인해 준다. 결합 에너지와 중력 질량의 비율은 약 0.24이며, 컴팩트도($\beta = M_g/R_{circ}$)와 선형 관계를 보인다.
  • 총 에너지: SQS의 총 에너지는 $\approx 5.5 \times 10^{47}$ J ($\approx 10^{54}$ ergs)로 추정된다. 총 에너지에 대한 외부 자기 에너지의 기여도는 1% 미만이다.
  • 안정성 기준: 축대칭 섭동에 대한 안정성은 $(dM/dR)_J > 0$을 보장함으로써 평가되었다.

• 관측과의 비교:

  • 모델이 예측하는 질량은 관측된 컴팩트 천체들과 호환된다. 특히, 비회전, 비자기 모델의 질량(2.35 $M_\odot$)은 PSR J0952-0607과 일치한다.
  • 회전/자기 모델의 높은 질량 한계(최대 2.80 $M_\odot$)는 GW190814 사건의 이차 구성 성분의 질량 범위(2.5–2.67 $M_\odot$)와 일치한다.
  • 또한 본 모델은 마그네타 SGR J1745–2900 및 PSR J1614-2230, PSR J0348+0432와 같은 다른 거대 펄사들의 질량 및 반지름 추정치를 수용한다.

의의 및 주장
본 논문은 란다우 양자화가 적용된 밀도 의존적 MIT 백 모델을 사용하는 제안된 SQS 모델이, 표준 쿼크 물질의 틀을 벗어난 이례적인 물리학 없이도 최근의 천문학적 사건에서 관측된 높은 중력 질량을 성공적으로 재현한다고 주장한다. 저자들은 $10^{19}$ G 미만의 자기장은 (Chatterjee 등의 선행 연구에서 언급된 바와 같이) EOS의 강성(stiffness)을 크게 변화시키지는 않지만, 회전과 결합될 경우 항성의 구조적 파라미터(질량 및 형태)를 유의미하게 변화시킨다는 점을 강조한다.

본 연구는 빠르게 회전하고 강하게 자기화된 SQS가 최대 2.8 $M_\odot$의 질량을 지지할 수 있음을 입증하며, 이는 LIGO/Virgo가 탐지한 "질량 간극(mass gap)" 천체 및 무거운 펄사에 대한 잠재적 설명을 제공한다. 저자들은 최소 중력 질량 한계와 케플러 구성에 대한 탐구는 동반 논문에서 다루어질 것임을 명시하였으며, 강하게 자기화되고 빠르게 회전하는 별의 임계 구성을 정확하게 계산하는 것이 향후 정밀화를 위한 한계점임을 인정하였다.