[논문 리뷰] Physical and chemical vertical structure of magnetostatic accretion disks of young stars

S. A. Khaibrakhmanov, A. E. Dudorov외 2명|arXiv (Cornell University)|2022. 10. 23.

Astrophysics and Star Formation Studies참고 문헌 19인용 수 2

한 줄 요약

이 연구는 젊은 T 타우리 별 주위의 자기적 평형 상태에 있는 원형 원반의 수직 물리적 및 화학적 구조를 모델링하며, 항성 중력, 난류 가열, 복사 가열, 고착된 대규모 자기장의 영향을 포함한다. 자기압력 기울기가 수직 평형 상태를 초월해 원반 두께를 증가시킬 수 있음을 발견하였고, CN 분자의 농도가 광학적 표면 근처와 원반 대기에서 최대가 되며, 이곳에서 자기장 강도가 약 0.01 G에 이르는 것으로 나타났다. 이는 이러한 영역에서 CN 선의 제이먼 분리 현상을 통해 자기장을 측정할 수 있음을 시사한다.

ABSTRACT

The vertical structure of the accretion disks of young stars with fossil large-scale magnetic field is studied. The equations of magnetostatic equilibrium of the disk are solved taking into account the stellar gravity, gas and magnetic pressure, turbulent heating, and heating by stellar radiation. The modelled physical structure of the disk is used to simulate its chemical structure, in particular, to study the spatial distribution of CN molecules. The disk of the typical T Tauri star is considered. Simulations show that the temperature within the disk in the region $r<50$ au decreases with height and density profiles are steeper than in the isothermal case. Outside the `dead' zone, vertical profiles of the azimuthal component of the magnetic field are nonmonotonic, and the magnetic field strength maximum is reached within the disk. The magnetic pressure gradient can cause an increase in the disk thickness in comparison with the hydrostatic one. The CN molecule concentration is maximum near the photosphere and in the disk atmosphere where the magnetic field strength at chosen parameters is $\sim 0.01$ G. Measurements of the Zeeman splitting of CN lines in the submm range can be used to determine the magnetic field strength in these regions of accretion disks.

연구 동기 및 목표

• 고착된 대규모 자기장이 젊은 별 주위의 원반 수직 구조에 미치는 영향을 조사하기 위해.
• 항성 중력, 난류 가열, 복사 가열, 자기압력의 복합적 영향을 원반 구조에 모델링하기 위해.
• 특히 자기장 강도와의 관계에서 원반 내 CN 분자의 화학적 분포를 시뮬레이션하기 위해.
• submillimeter 범위의 CN 선 제이먼 분리를 통해 원반 광학적 표면 및 대기에서 자기장 강도를 측정할 수 있는지 타당성을 평가하기 위해.

제안 방법

• 원통좌표계(r, φ, z)에서 수직 중력, 기체 압력, 자기압력, 난류 가열을 포함한 자기적 평형 상태 방정식을 해결한다.
• 원반 표면에서 복사 가열을 위해 복사 플럭스 경계 조건을 적용하며, 흡수 분율 f = 0.05로 고정한다.
• 난류 점성에 대해 Shakura-Sunyaev의 α-예측법을 사용하고, 유도 방정식에 옴 손실, 이온-중성자 분리, 홀 효과를 포함한다.
• 다섯 개의 경계 조건을 사용하여 방정식 시스템을 해결한다: z=0에서 대칭 조건, z=z_s에서 효과적 온도, 원반 중면에서 플럭스 및 자세각 자기장이 0.
• DVODE 적분기로 10^6년 동안 화학적 진화 모델링을 수행하며, 화학적 균형 방정식 시스템을 사용한다.
• r = 0.25, 10, 및 50 au에서 밀도, 온도, 자기장(Bφ), 그리고 수소에 대한 CN 농도(x_CN)의 수직 프로파일을 계산한다.

실험 결과

연구 질문

• RQ1 고착된 대규모 자기장 존재 시, 절연체 또는 수직 평형 상태에 비해 원형 원반의 수직 온도 및 밀도 구조는 어떻게 변화하는가?
• RQ2 원반 내에서 자세각 자기장 성분(Bφ)의 공간 분포는 어떻게 되며, 높이에 따라 비단조화적인 행동을 보이는가?
• RQ3 원반 내에서 CN 분자의 농도가 어디에서 최대가 되는가? 그리고 어떤 물리적 조건이 그 형성 또는 생존을 촉진하는가?
• RQ4 submillimeter 범위의 CN 선 제이먼 분리를 통해 원반 광학적 표면 및 대기에서 자기장 강도를 추정할 수 있는가?

주요 결과

• 내부 원반(r < 50 au)에서는 온도가 높이에 따라 감소하여, 등온 상태에 비해 더 급격한 밀도 프로파일을 보이며, 광학적 표면 근처의 밀도는 수직 평형 상태에 비해 4–5개의 주자 정도 감소한다.
• 자세각 자기장(Bφ)은 원반 내부에서 최대값을 보이며(예: r = 50 au일 때 z ≈ 2H에서 약 0.01 G), '사망 영역' 외부에서는 비단조화적인 수직 프로파일을 나타낸다.
• 자기압력 기울기가 수직 평형 상태를 초월해 효과적인 원반 두께를 증가시킬 수 있으며, 특히 Bφ가 내부에서 최대가 되는 영역에서 두드러진다.
• CN 분자의 농도(x_CN)는 광학적 표면(z ≈ z_s) 근처와 상부 대기에서 최대값인 약 10⁻¹⁰에서 10⁻⁹에 이르며, 이는 이전 관측 결과와 일치한다.
• 이 높은 CN 농도 영역의 자기장 강도는 약 0.01 G이며, 이는 submillimeter 범위의 CN 선 제이먼 분리를 통해 직접 자기장을 측정할 수 있음을 의미한다.
• 모델은 항성 복사에 의한 원반 가열로 광학적 표면에서 급격한 온도 상승이 발생하며, 이로 인해 상부 층에서 난류 가열에서 복사 가열으로의 전이가 일어남을 예측한다.