A Three-Dimensional Exploration of Magnetic Fields, Rotation, and Shock Revival in a $39 M_\odot$ Core-Collapse Supernova Progenitor

기술적 요약: 39 M⊙ 핵붕괴 초신성 progenitor 에서의 자기장, 회전, 충격파 부활에 대한 3 차원 탐구

문제 및 동기
핵붕괴 초신성 (CCSNe) 에서 정체된 충격파의 부활은 중성미자 가열, 회전, 자기장, 그리고 유체역학의 상호작용을 포함하는 복잡한 문제로 남아 있습니다. 중성미자 메커니즘은 특히 높은 컴팩트성을 가진 많은 progenitor 에서 견고하지만, 자기회전 (MR) 메커니즘은 충격파 부활과 고에너지 천체 현상 (예: 초신성, Type Ic-BL 초신성) 을 구동하는 대안적 경로를 제공합니다. 이 메커니즘은 와동 (winding) 과 불안정성을 통해 회전 자유 에너지를 자기 에너지로 변환하고 쌍극성 유출을 구동하기 위해 빠른 회전과 강력한 자기장을 필요로 합니다. 그러나 이전 연구들은 종종 축대칭 (2D) 에 의존해 왔는데, 이는 제트를 인위적으로 안정화시키고 키크 불안정성과 같은 비축대칭 불안정성을 억제할 수 있습니다. 더 나아가, 높은 컴팩트성과 빠른 회전을 가진 3 차원 환경에서 자기장이 충격파 부활 시간 척도, 유출 형태, 그리고 블랙홀 형성에 대한 원시 중성자별 (PNS) 의 안정성에 미치는 구체적인 영향은 추가적인 조사가 필요합니다.

방법론
저자들은 높은 컴팩트성 ($\xi_{2.5} = 0.75$) 과 빠른 회전 ($\Omega_c \approx 0.58$ rad s$^{-1}$) 을 가진 39 M⊙ 제로 주계열 (ZAMS) progenitor 의 3 차원 유체역학 (HD) 및 자기유체역학 (MHD) 시뮬레이션을 수행했습니다. Progenitor 는 MESA 를 사용하여 준화학적으로 균질한 진화로 진화시켰습니다. 세 가지 서로 다른 모델을 FLASH 코드와 Spark 솔버를 사용하여 시뮬레이션했으며, 적응형 메쉬 정밀화 (AMR) 와 M1 닫힘을 갖춘 다차원, 에너지 의존적, 2-모멘트 중성미자 수송 방식을 사용했습니다:

1. B1R1: 회전하는, 자기화된 모델. 회전 축과 정렬된 극성 시드 자기장 ($B_0 \approx 5 \times 10^{11}$ G) 을 초기화했습니다.
2. B0R1: 회전하는, 비자기화된 제어 모델 (HD 만).
3. B0R0: 비회전, 비자기화된 기준 모델 (HD 만).

시뮬레이션은 붕괴 전부터 코어 반동 (core bounce) 을 거쳐 초기 반동 후 단계 (약 1 초까지) 로 진화했으며, 코어 근처에서 최상 해상도 약 0.61 km 를 가진 직교 격자를 사용했습니다. 일반 상대론적 효과는 Marek 등 (2006) 의 Case A 를 사용한 유효 퍼텐셜로 근사화되었습니다.

주요 결과
세 모델 모두 시뮬레이션 구간 내에서 runaway 충격파 확장을 달성했지만, 형태, 시간 척도, 그리고 잔해의 운명에서 현저하게 다른 특징을 보였습니다:

• 충격파 부활 및 형태:

  • B1R1 (자기화/회전): 가장 빠른 충격파 부활 ($t_{rev} \approx 150$ ms, 반동 후) 을 보였으며 가장 명확한 쌍극성 유출을 발달시켰습니다. 유출은 초기에 제트와 유사했으나 비축대칭 불안정성에 의해 확장되고 왜곡되었습니다.
  • B0R0 (비회전): 강한 중성미자 구동 난류를 통해 더 늦게 ($t_{rev} \approx 250$ ms) 부활했습니다. 충격파는 선호되는 극 방향 없이 큰 무질서한 플룸으로 확장되었습니다. 중요하게도, 이 모델은 회전 지지의 부재로 인해 반동 후 약 1 초 후에 블랙홀로 붕괴했습니다.
  • B0R1 (회전/비자기화): 가장 늦은 충격파 부활 ($t_{rev} \approx 370$ ms) 을 경험했습니다. 회전은 대류를 억제하고 중성미자 광도를 감소시켜 B0R0 에 비해 더 컴팩트한 진화를 이끌었습니다. 명확한 쌍극성 유출은 발달하지 않았습니다.

• 각운동량 및 에너지 수송:

  • B1R1에서 맥스웰 응력은 내부 PNS 에서 바깥쪽으로 각운동량을 효율적으로 재분배하여 내부 코어의 스핀을 감소시키고 회전 자유 에너지를 추출했습니다. 이 과정은 차등 회전의 저장소를 활용하여 극 유출을 구동했습니다.
  • 각운동량 예산은 대류가 내부 코어에 각운동량을 공급했지만, 자기 토크가 스핀 업을 상쇄하는 지배적인 싱크 역할을 했음을 보여주었습니다. 중성미자 방출은 2 차적이고 더 작은 싱크로 작용했습니다.
  • B1R1 의 진단 에너지 ($E_{diag}$) 는 회전 에너지 추출과 일치하여 가장 일찍 그리고 가장 빠르게 증가했습니다. 자기 에너지 예산은 추출된 에너지의 상당 부분이 해결된 자기 에너지로 유지되는 대신, 미해결 재결합 및 서브-그리드 전이와 관련된 유효 소산 항 ($\dot{Q}_{diss}$) 을 통과했음을 나타냈습니다.

• 불안정성 및 3 차원 역학:

  • B1R1 의 자기화 유출은 비축대칭 키크 불안정성 ($m=1$ 모드) 의 영향을 받았습니다. 진단 결과 자기 에너지 무게중심 변위의 초기 지수적 성장 (e-곱 시간 약 1.5 ms) 을 보였으며, 이는 이론적 키크 시간 척도와 일치합니다.
  • 안정적이고 좁은 제트를 예측할 수 있는 축대칭 시뮬레이션과 달리, 3 차원 시뮬레이션은 유출 깔때기가 시간이 지남에 따라 기울어지고 (몇 도) 확장됨을 보여주었습니다. 유출은 생존하여 큰 반경까지 전파되었으나 좁은 콜리메이션을 잃었으며, 이는 제트 안정성을 정확하게 모델링하기 위해 3 차원 효과가 중요함을 시사합니다.

• 핵합성:

  • 분출물의 사후 처리는 뚜렷한 조성 차이를 보여주었습니다. B0R0은 중성미자 구동 폭발에 전형적인 양성자 과잉 분출물 ($Y_e > 0.5$) 을 생성하여 철군 원소 ($A \sim 56$) 까지 합성했습니다.
  • B1R1은 빠른 유출 발달로 인해 약한 r-과정을 가능하게 하는 상당한 중성자 과잉 성분 ($Y_e$ 약 0.33 까지) 을 생성하여 $A \sim 90$까지 합성했습니다. 그러나 초기 시드 자기장과 키크에 의한 확장은 더 극단적인 자기장을 가진 모델에서 보이는 극도로 중성자 과잉 분출물 ($Y_e < 0.32$) 의 생성을 방지했습니다.

의의 및 주장
이 논문은 극단적이고 높은 컴팩트성을 가진 progenitor 에서 빠른 회전과 강력한 자기장의 결합이 본질적으로 3 차원인 초기 자기장 보조 충격파 부활 경로를 열 수 있다고 주장합니다. 주요 결론은 다음과 같습니다:

1. 구동자로서의 자기장: 자기 응력은 단순히 유출을 형성하는 것이 아니라, 회전 자유 에너지를 적극적으로 추출하고 내부 코어의 스핀을 감소시키며, 비자기화 회전 경우보다 충격파 부활을 촉진합니다.
2. 3 차원의 필요성: 자기장 구동 유출의 진화는 본질적으로 비축대칭적입니다. 키크 불안정성은 유출을 확장하고 기울게 하여 지속적이고 좁은 제트의 형성을 방지합니다. 축대칭 시뮬레이션은 이러한 유출의 안정성과 콜리메이션을 과대평가할 수 있습니다.
3. 잔해의 안정성: 회전은 높은 컴팩트성 progenitor 에서 즉각적인 블랙홀 형성에 대한 상당한 지지를 제공하지만, 시뮬레이션 구간을 넘어선 잔해의 장기적 안정성은 여전히 불확실합니다.
4. 에너지 경로: 차등 회전에서 유출로의 에너지 전달은 상당한 미해결 소산을 수반하며, 이는 MHD CCSN 시뮬레이션에서 서브-그리드 물리학의 중요성을 강조합니다.

저자들은 미해결 소규모 다이나모 증폭을 보상하기 위해 강한 시드 자기장을 사용했다는 점, 소산 진단에서의 해상도 민감성, 그리고 회전 모델의 궁극적인 운명을 결정하기 위한 더 긴 시뮬레이션의 필요성을 포함한 제한 사항을 지적했습니다. 그들은 시뮬레이션이 자기회전 충격파 부활의 그림을 지지하지만, 구체적인 결과는 유출의 전체 3 차원 역학에 크게 의존한다고 강조합니다.