[논문 리뷰] Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal MHD case

Jacques Masson, G. Chabrier외 3명|Figshare|2015. 09. 18.

Astrophysics and Star Formation Studies참고 문헌 51인용 수 79

한 줄 요약

이 논문은 저질량 항성 형성에서 비이상적 자기유체역학(MHD) 시뮬레이션을 사용하여 이중극자 확산의 영향을 조사하며, 이는 자기장 확산 장벽을 통해 자기장 농축을 억제하고, 이로 인해 킬러피안 원반의 형성이 가능해짐을 보여준다. 이는 이상적 MHD에서는 각운동량이 과도하게 억제되어 원반 형성이 불가능한 것과 대조된다. 이중극자 확산을 포함한 모델은 초기 자기장의 도핑 정도에 따라 질량 0.02–0.2 M⊙, 반경 20–80 au 범위의 물리적으로 타당한 원반을 형성한다.

ABSTRACT

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).

연구 동기 및 목표

• 이상적 MHD 시뮬레이션에서의 비물리적인 자기장 농축과 과도한 자기장 브레이킹 문제를 해결하기 위해.
• 이중극자 확산이 첫 번째 라슨 코어 형성 과정 중 각운동량 이동과 자기장 농축 재분포를 어떻게 조절하는지 조사하기 위해.
• 비이상적 MHD가 클래스-0 및 클래스-I 원반성 항성의 관측 결과와 일치하는 회전으로 지지되는 원반을 생성할 수 있는지 확인하기 위해.
• 초기 자기장 방향과 자기장 도핑 정도가 원반 형성과 안정성에 미치는 영향을 평가하기 위해.
• 수치적 방법이 자기확산을 시뮬레이션할 때 수치적 확산에 의한 비물리적 아티팩트를 피하는 데서의 역할을 평가하기 위해.

제안 방법

• 1 M⊙ 전성운핵의 붕괴를 시뮬레이션하기 위해 비이상적 MHD 해법을 갖춘 RAMSES 적응 격자 분해 코드를 사용하였다.
• 이온, 중성, 먼지 종류의 농도를 고려한 단순화된 화학 네트워크를 도입하여 확산 계수를 계산함으로써 이중극자 확산을 포함시켰다.
• 핵 근처에 자기장 확산 장벽(분리 단계)을 적용하여 자기장 강도가 약 0.1 G를 초과하지 않도록 하였다.
• 수치적 확산 영향을 시험하기 위해 다양한 리만 해법(HLL, HLLD, Lax-Friedrich)을 사용하였다.
• 원반 특성을 반경 방향 속도 프로파일을 기반으로 정의하고, 비킬러피안 편이를 보정한 후 킬러피안 기대치와 비교하였다.
• 초기 자기장 도핑 정도(μ = 2에서 μ = 5)와 자기장-회전축 간의 기울기 각도를 변화시켜 시뮬레이션을 수행하였다.

실험 결과

연구 질문

• RQ1 붕괴하는 전성운핵에서 이상적 MHD와 비교해 이중극자 확산은 자기장 재분포에 어떻게 영향을 미치는가?
• RQ2 비이상적 MHD 시뮬레이션은 클래스-0 및 클래스-I 원반성 항성의 관측 결과와 일치하는 회전으로 지지되는 원반을 생성할 수 있는가?
• RQ3 자기장 확산 장벽이 자기장 강도 증폭을 조절하고 과도한 브레이킹을 방지하는 데서 어떤 역할을 하는가?
• RQ4 수치적 해법 선택(HLL 대비 HLLD 등)은 첫 번째 라슨 코어의 진화와 원반 형성에 어떤 영향을 미치는가?
• RQ5 초기 자기장 방향과 자기장 도핑 정도가 비이상적 MHD 시뮬레이션에서 원반 질량과 크기에 얼마나 영향을 미치는가?

주요 결과

• 이중극자 확산은 첫 번째 라슨 코어 근처에 자기장 확산 장벽을 형성하여 자기장 강도 증폭을 약 0.1 G로 제한하고, 자기장 농축을 방지한다.
• 원반 질량은 초기 자기장 도핑 정도에 따라 약 0.02 M⊙(μ = 2)에서 0.2 M⊙(μ = 5) 사이로 변동하며, 반경은 20–30 au에서 40–80 au 사이로 변화한다.
• 모든 비이상적 MHD 시뮬레이션에서 킬러피안 속도 프로파일이 회복되며, 이는 이상적 MHD에서는 원반 형성이 억제됨을 보여준다.
• 이상적 MHD에서는 강력한 자기장 외상류가 관찰되나, 비이상적 MHD에서는 자기장 브레이킹 감소로 인해 약화되거나 사라진다.
• 비이상적 MHD에서는 자기장과 회전축 간의 기울기가 원반 형성에 거의 영향을 주지 않으며, 이는 이상적 MHD와는 대조적으로 물리적 영향이 우세함을 확인한다.
• 이상적 MHD 시뮬레이션은 반대 방향 회전 외상류, 교환 불안정성, 수치적 자기장 재분포 등 비물리적 행동을 보이며, 비이상적 MHD에서는 이러한 현상이 관찰되지 않는다.