[논문 리뷰] The role of turbulence and winding in the development of large-scale, strong magnetic fields in long-lived remnants of binary neutron star mergers

Ricard Aguilera-Miret, Carlos Palenzuela외 2명|arXiv (Cornell University)|2023. 07. 10.

Gamma-ray bursts and supernovae인용 수 4

한 줄 요약

이 연구는 고해상도 대류모의(LES)와 기울기 미소스케일 모델을 사용하여 장수명 이중 중성자별 병합 잔여물의 자기장 진화를 조사한다. 결과적으로 약 30ms 후에 자기장 감도에 의해 토로이달 성분이 선형적으로 증가하는 것을 발견하였으며, 다항 성분은 약 90ms 이후에야 소규모에서 성장하기 시작한다. 이는 시뮬레이션 시간스케일 내에서 대규모 다항 성분 또는 제트가 형성되지 않음을 시사하며, 초기 강력한 다항 자기장은 현실적인 병합 후 구조와 부합하지 않는다.

ABSTRACT

We perform a long and accurate Large-Eddy Simulation of a binary neutron star merger, following the newly formed remnant up to 110 milliseconds. The combination of high-order schemes, high-resolution and the gradient subgrid-scale model allow us to have among the highest effective resolutions ever achieved. Our results show that, although the magnetic fields are strongly amplified by the Kelvin-Helmholtz instability, they are coherent only over very short spatial scales until t \gtrsim 30 ms. Around that time, magnetic winding becomes more efficient leading to a linear growth of the toroidal component and slowly ordering the field to more axisymmetric, large scales. The poloidal component only starts to grow at small scales at much later times t \gtrsim 90 ms, in a way compatible with the magneto-rotational instability. No strong large-scale poloidal field or jet is produced in the timescales spanned by our simulation, although there is an helicoidal structure gradually developing at late times. We highlight that soon after the merger the topology is always strongly dominated by toroidal structures, with a complex distribution in the meridional plane and highly turbulent perturbations. Thus, starting with strong purely dipolar fields before the merger is largely inconsistent with the outcomes of a realistic evolution. Finally, we confirm the universality of the evolved topology, even when starting with very different magnetic fields confined to the outermost layers of each neutron star.

연구 동기 및 목표

• 장수명 이중 중성자별 병합 잔여물에서 대규모 강력한 자기장의 발달을 조사하기 위해.
• 병합 후 자기장 구조를 형성하는 난류와 자기장 감도의 역할을 평가하기 위해.
• 초기 강력한 대규모 다항 자기장이 현실적인 병합 동역학과 물리적으로 일치하는지 평가하기 위해.
• 다양한 초기 자기장 구성과 수치 해상도에서 자기장 진화의 탄력성 테스트를 위해.
• 고해상도 LES와 기울기 SGS 모델링이 소규모 난류가 해결되지 않은 상태에서도 핵심 MHD 불안정성과 자기장 진화를 포착할 수 있는지 평가하기 위해.

제안 방법

• 110ms까지의 장기적 정밀 대류모의(LES)를 통해 이중 중성자별 병합을 수행하기 위해.
• 고차수 수치 스킴과 기울기 미소스케일(SGS) 모델을 사용하여 효과적 해상도를 향상시키고 해결되지 않은 MHD 역학을 포착하기 위해.
• 현실적인 상태방정식을 사용하여 병합를 시뮬레이션하고, 초기 구성에서 자기장을 외부 층에 국한하거나 별 전체에 분포시켜 시뮬레이션하기 위해.
• 다항 및 토로이달 자기장 성분의 진화를 추적하며, 응집도, 성장률, 구조적 특성에 중점을 두기 위해.
• 켈빈-헬름홀츠 불안정성(KHI), 자기장 감도, 자기장-회전 불안정성(MRI) 간의 상호작용 분석하기 위해.
• 다양한 초기 자기장 구성과 수치 해상도 간의 결과 비교를 통해 수렴성과 탄력성 평가하기 위해.

Figure 1: Magnetic field evolution in the orbital plane . The magnetic field strength $|\vec{B}|$ is shown in logarithmic scale at times $t=\{0.5,1.25,2.5,3.75,5,15,25,50,110\}$ ms after the merger. Outer and inner black lines depict density contours of $\rho=10^{13}~{}\rm{g~{}cm^{-3}}$ and $5\times

실험 결과

연구 질문

• RQ1 난류와 자기장 감도는 장수명 이중 중성자별 병합 잔여물에서 대규모 자기장의 성장과 응집도에 어떻게 영향을 미치는가?
• RQ2 초기 자기장 구성(예: 다항 대비 표면에 국한된)이 최종 자기장 구조와 진화에 어느 정도의 영향을 미치는가?
• RQ3 LES에서 기울기 SGS 모델이 KHI와 MRI와 같은 해결되지 않은 MHD 불안정성의 본질적 역학을 정확히 포착할 수 있는가?
• RQ4 강력한 초기 다항 자기장이 포함된 시뮬레이션에서 왜 현실적인 병합 잔여물에서 관찰되는 토로이달 우세성과 일치하지 않는가?
• RQ5 대규모 다항 자기장은 어느 시간스케일에서 성장하기 시작하며, 이는 제트 형성 메커니즘과 호환되는가?

주요 결과

• 켈빈-헬름홀츠 불안정성에 의해 자기장은 강하게 증폭되지만, 응집도는 약 30ms 이상이 되기 전까지 몇 킬로미터 규모에 국한된다.
• 약 30ms 이후로 자기장 감도가 지배적이 되어 토로이달 자기장 성분이 시간에 비례하여 선형적으로 증가하고 축대칭적이고 대규모의 정렬을 촉진한다.
• 다항 자기장 성분은 약 90ms 이전까지 강도가 거의 일정하게 유지되며, 이후 소규모에서 성장하기 시작하여 자기장-회전 불안정성(MRI)과 일치한다.
• 110ms 이내의 시뮬레이션 시간스케일 내에서는 강력한 대규모 다항 자기장이나 상대론적 제트가 형성되지 않았지만, 점차 나선형 구조가 발달된다.
• 초기 순수한 다항 자기장은 병합 후 구조와 부합하지 않으며, 어떤 초기 구성이든 병합 후에는 난류적이고 토로이달적인 구조가 지배한다.
• 자기장 증폭과 진화는 초기 자기장 분포 방식(외부 층에 국한된 경우 또는 별 전체에 분포된 경우)에 관계없이 강건하게 유지되며, 이는 보편적인 자기장 구조 진화를 시사한다.

Figure 2: 2D plots in the meridional plane . The rows show, from top to bottom, the rest-mass density in [g cm -3 ], magnetic field intensity in [G], angular velocity of the fluid $\Omega$ in [rad s -1 ], and inverse of $\beta$ factor, at $t=\{10,50,100\}$ ms after the merger (first, second and thir