$\Delta(1232)$ isobars and weak-interaction bulk viscosity in neutron-star merger matter

기술 요약: $\Delta(1232)$ 아이소바와 중성자별 병합 물질 내 약한 상호작용 벌크 점성

문제 제기
이중 중성자별 병합은 핵 밀도의 몇 배($n_0$)에 달하는 밀도와 수십 MeV의 온도를 가진 핵 물질을 탐사한다. 이러한 조건 하에서는 더 무거운 바리온인 $\Delta(1232)$ 아이소바가 나타날 수 있다. 기존의 병합 물질 내 벌크 점성 계산은 핵자 물질(뮤온 포함 또는 미포함)과 하이퍼론 물질을 고려해 왔으나, $\Delta$ 아이소바의 구체적인 역할은 다루어지지 않았다. 강한 상호작용($NN \leftrightarrow N\Delta$)에 의해 강한 평형 상태를 유지하는 $\Delta$ 아이소바의 존재가 맛(flavor) 평형을 담당하는 약한 상호작용율과 그로 인한 밀도 진동의 벌크 점성 감쇠를 어떻게 변화시키는지에 대한 이해에 결정적인 공백이 존재한다.

방법론
저자는 강한 평형 상태에 있는 $\Delta$ 쿼텟($\Delta^{--}, \Delta^-, \Delta^0, \Delta^+$)을 포함하는 뉴트리노 투과성 물질에 대해 기존의 완전 위상 공간 계산을 확장하였다. 본 연구는 두 가지 상대론적 평균장 모델(QMC-RMF3 및 DDME2)을 사용하며, 이 모델들은 차가운 $\beta$-평형 상태의 핵자 직접 우르카 과정(direct Urca process)을 $8n_0$까지까지 금지하도록 설계되어 있어, 이 과정이 열리는 것은 오직 $\Delta$ 아이소바에 의한 것임을 보장한다.

주요 방법론적 구성 요소는 다음과 같다:

• 행렬 요소(Matrix Elements): Rarita–Schwinger 스핀 합을 사용하여 $\Delta \to N$ 및 $\Delta \to \Delta'$의 전하 변화 전이에 대한 상대론적, 스핀 합산 행렬 요소를 유도한다. 본 연구는 모든 전자 및 뮤온 우르카 채널을 포함한다.
• 율 계산(Rate Calculation): NWA(Nucleon Width Approximation)를 사용하여 전체 위상 공간에서 우르카 과정을 평가한다. 이 접근법은 충돌 확산("수정된 우르카")을 직접 우르카 기여와 동일한 방식으로 취급한다.
• 매질 내 폭(In-Medium Widths): 채워진 $\Delta$ 상태가 강한 붕괴에 대해 갖는 안정성을 결정하기 위해 파울리 차단(Pauli blocking)을 고려한 매질 내 $\Delta \to N\pi$ 폭을 계산한다.
• 선형 응답(Linear Response): 유한 온도 맛 평형에 대한 전자 및 뮤온 분율($Y_e, Y_\mu$)의 결합된 선형 응답을 통해 벌크 점성을 공식화한다. 완화 행렬(relaxation matrix)은 서로 다른 전하를 가진 상태로 바리온을 변환하는 모든 과정을 포함한다.
• 검증: 기존 문헌의 핵자 결과를 재현하고, 열역학적 항등식(Gibbs–Duhem 관계식) 및 행렬 요소 합을 명시적인 스피너 계산과 대조하여 코드를 검증하였다.

주요 기여 및 결과

1. 조성 변화에 의한 직접 우르카 개방:
   주요 발견은 $\Delta^-$ 아이소바가 임계치(모델에 따라 $1.45n_0$에서 $2.3n_0$ 사이)를 넘어서면 전하 중성을 유지하는 데 있어 렙톤의 역할을 대신한다는 것이다. 이러한 대체는 전자 분율을 10배 가량 낮추고 양성자 분율을 실질적으로 두 배 증가시킨다. 결과적으로, 핵자 직접 우르카 과정($n \to p e \bar{\nu}_e$)의 운동학적 조건이 충족되어, 선택된 핵자 모델에서 이전에 닫혀 있던 채널이 열리게 된다.

2. 2단계 평형 메커니즘:
   본 연구는 두 가지 뚜렷한 평형 영역을 식별한다:

• 저온 ($\lesssim 3$ MeV): 조성 변화(핵자 직접 우르카 과정의 개방)가 평형 속도를 제어한다. 이 영역에서 $\Delta$ 약한 과정 자체는 운동학적으로 금지되어 있으며 기여도가 미미하다.
• 고온 ($\gtrsim 3.5$ MeV): 열적 흐림(thermal blurring)과 충돌 확산으로 인해 $\Delta$ 약한 과정($\Delta \to N l \bar{\nu}$ 및 $\Delta \to \Delta' l \bar{\nu}$)이 활성화된다. 이 영역에서 $\Delta$ 채널은 전체 평형 속도의 절반 이상을 차지한다.

3. 벌크 점성 공명 이동:
   $\Delta$ 아이소바의 존재는 벌크 점성($\zeta$)을 크게 변화시킨다. $\Delta$ 유도 조성 변화와 그에 따른 약한 과정이 완화 속도를 가속화하기 때문에, 벌크 점성의 공명 극대값(완화 속도가 진동 주파수와 일치하는 지점)이 더 낮은 온도로 이동한다.

• $1$kHz 진동의 경우, 순수 핵자 물질에서의 공명 온도는 약$4\text{--}6$MeV에서$\Delta$혼합 물질에서는 약$2\text{--}3.7$ MeV로 이동한다.
• 공명 미만(Below Resonance): $2\text{--}3$ MeV의 온도에서 완화 속도의 가속은 시스템을 공명에 가깝게 만들어, 핵자 물질에 비해 벌크 점성을 10배 이상 증가시킨다.
• 공명 초과(Above Resonance): $5\text{--}10$ MeV 및 고밀도에서의 높은 온도에서는 시스템이 공명보다 훨씬 높은 상태로 밀려나, 벌크 점성이 1~3 자릿수만큼 감소한다.

4. kHz 진동의 감쇠:
   밀도 진동이 $20$ ms 이내(병합 초기 단계와 관련된 시간 척도)에 감쇠되는 영역이 크게 확장된다. 핵자 물질에서 $1$kHz 진동에 대한 효율적인 감쇠는 포화 밀도 근처로 제한된다.$\Delta$아이소바가 있으면 효율적인 감쇠 영역은 포화 밀도에서$3\text{--}5n_0$까지, 그리고$2\text{--}2.75$ MeV까지 낮아진다. 가장 짧은 감쇠 시간은 핵자 물질보다 몇 배 더 짧다.

5. 평형 상태에 대한 민감도:
   본 논문은 율(rate)을 선형화하는 기준이 되는 평형 상태의 선택이 매우 중요함을 보여준다. 실제 유한 온도 평형이 아닌 차가운 $\beta$-평형 조건(화학적 불균형이 0인 상태)에 대해 선형화하면 고밀도에서 공명 온도가 최대 몇 MeV까지 이동한다. 이는 차가운 근사를 사용했던 이전 연구들과의 불일치를 설명해 준다.

6. $\Delta$ 아이소바의 안정성:
   $\Delta \to N\pi$ 붕괴가 운동학적으로 허용됨에도 불구하고, 채워진 $\Delta$ 아이소바는 밀도가 높은 매질 내에서 강한 붕괴에 대해 파울리 차단된다. 이들은 강한 붕괴보다는 충돌에 의해 결정되는 스펙트럼 폭을 가진 안정적인 준입자(quasiparticle)처럼 행동한다.

의의 및 주장
본 논문은 중성자별 병합 물질 내 $\Delta$ 아이소바의 존재가 주로 조성 변화를 통해 핵자 직접 우르카 과정을 개방함으로써, 밀도 진동의 벌크 점성 감쇠를 근본적으로 변화시킨다고 주장한다. 이 효과는 수 MeV의 온도에서 약한 상호작용 응답을 지배하며, 이는 병합 후 단계와 관련된 온도이다.

저자는 다음과 같은 사실을 시사한다:

• 더 낮은 온도 영역에서의 감쇠 강화: 만약 $\Delta$ 아이소바가 존재한다면, 핵자 물질로는 진동을 유의미하게 감쇠하기에 너무 느린 더 낮은 온도 영역($2\text{--}3$ MeV)에서도 효과적인 벌크 점성 감쇠가 발생할 수 있다.
• 모델 의존성: 이러한 효과의 크기는 특정 상태 방정식과 $\Delta$ 결합 매개변수에 따라 달라지지만, 공명 온도의 질적인 이동은 테스트된 모델 전반에서 견고하게 나타난다.
• 방법론적 필요성: 정확한 병합 역학 모델링을 위해서는 약한 상호작용 율을 차가운 극한이 아닌 유한 온도 평형에 대해 선형화해야 하며, 이 차이는 예측된 공명 온도에 상당한 영향을 미친다.

본 연구는 $\Delta$ 아이소바가 새로운 느린 완화 변수를 도입하지는 않지만(강한 평형 상태이므로), 열역학적 응답과 약한 과정의 집합을 모두 수정하여 벌크 점성을 위한 더 빠른 평형과 더 낮은 온도에서의 공명 이동을 초래한다고 결론짓는다.