AB Aur, a Rosetta stone for studies of planet formation (IV): C/O estimates from CS and SO interferometric observations

기술 요약: AB Aur, 행성 형성 연구의 로제타 스톤 (IV): CS 및 SO 간섭계 관측을 통한 C/O 추정

문제 및 배경
원시행성계 원반은 행성의 탄생지이며, 행성의 대기 조성을 결정하는 초기 화학적 함량을 설정한다. 이 맥락에서 중요한 매개변수는 탄소-산소(C/O) 비율로, 이는 행성의 광물학 및 대기 화학에 영향을 미친다. 황(Sulfur)은 우주에서 열 번째로 풍부한 원소이지만, 성간 매질(ISM) 및 원시행성계 원반에서의 화학적 거동은 여전히 잘 이해되지 않았으며, 밀도가 높은 가스 내에서 우주 표준값에 비해 현저히 결핍된 모습을 보이곤 한다. 또한, 행성은 원반의 중간면(mid-plane)에서 형성되지만, 대부분의 분자 추적자(tracer)는 상층부의 "따뜻한 분자층(warm molecular layer)"에서 기원하므로, 행성 형성 영역을 연구하는 데 어려움이 따른다. 본 연구는 이전에는 SO 검출의 희소성으로 인해 제한적이었던 CS/SO 함량 비를 활용하여, AB Aur 원시행성계 원반의 황 함량과 C/O 비율을 규명하고자 한다.

방법론
저자들은 2022년 초 IRAM NOEMA 어레이(Project W21BA)를 통해 얻은 AB Aur 방향의 CS 3-2 전이에 대한 새로운 간섭계 관측 데이터를 제시한다. 이 데이터는 SO(Project 2021.1.00690)에 대한 기존 ALMA 관측 및 다른 분자 종(HCN, HCO+, CO 동위원소, H2CO, H2S)에 대한 이전의 NOEMA 관측 데이터와 결합되었다.

분석은 엄격한 다단계 과정을 따랐다:

1. 운동학적 특성 규명: 속도 적분 강도 지도(moment 0), 1차 모멘트 지도(강도 가중 속도), 2차 모멘트 지도를 생성하였다. inclination($i$)과 position angle(PA)은 eddy 코드(Teague 2019)와 bettermoments를 사용하여 여러 종의 1차 모멘트 지트를 피팅함으로써 도출되었다.
2. 케플러 마스킹(Keplerian Masking): 도출된 원반 파라미터($i \approx 22^\circ$, $PA \approx 237^\circ$, $M_* = 2.4 M_\odot$, $d = 155.9$ pc)를 사용하여 케플러 운동을 하는 방출을 비케플러 운동(예: 나선형 구조, 스트리머)으로부터 분리하기 위한 케플러 마스크를 구축하였다. 이를 통해 깨끗한 칼럼 밀도 지도를 생성하고 잔여 비케플러 특징들을 식별할 수 있었다.
3. 칼럼 밀도 추정: 국소 열역학적 평형(LTE)과 가스 온도 30 K에서의 광학적으로 얇은 방출을 가정하여, CS와 SO의 칼럼 밀도를 계산하였다.
4. 화학 모델링: 관측 데이터를 세 가지 핵심 파라미터(가스-먼지 질량비, 황 함량 [S/H], C/O 비율)를 변화시킨 NAUTILUS 화학 모델(Rivière-Marichalar et al. 2020) 세트와 비교하였다. 모델은 공동(cavity, $R_{in}=100$ au, $R_{out}=300$ au)을 포함하도록 조정되었으며, 관측된 빔 크기에 맞춰 컨볼루션(convolution)되었다.

주요 결과

• 원반 기하구조: 분석 결과, 모든 관측된 분자 종에서 일관되게 나타나는 원반 경사각 $22.0^\circ \pm 0.5^\circ$ 및 위치각 $237.0^\circ \pm 0.7^\circ$를 확인하였다.
• 공간 분포: CS 3-2 방출은 SO($\sim 1.4''$)와 H2S($\sim 1.2''$)의 피크 사이에 위치한 $\sim 1.3''$ ($\sim 212$ au)에서 피크를 이루는 고리 형태를 띤다. CS 고리의 폭(FWHM)은 약 $0.86''$ ($\sim 141$ au)이다.
• 운동학: CS, H2S, SO는 지배적으로 케플러 운동을 보이며 유의미한 비케플러 잔차를 보이지 않았다($3\sigma$ 수준). 반면, CO 동위원소는 대규모 나선 팔과 스트리머를 추적하는 강한 비케플러 잔차를 보였고, HCO+는 정렬이 어긋난 내부 원반 및 후보 행성인 AB Aur b와 관련된 잔차를 보였다. 이는 CS와 H2S가 CO보다 정적인 중간면을 더 잘 추적하는 지표임을 시사한다.
• 화학적 분리: 조사된 종들의 방사형 방출 피크와 결합 에너지 사이의 상관관계가 발견되지 않았으며, 이는 열적 탈착(thermal desorption)이 이 원반의 화학적 분리를 일으키는 주요 동인이 아님을 나타낸다.
• C/O 비율 및 황 함량:
  • 관측된 CS/SO 칼럼 밀도 비는 1.8에서 2.6 사이이며, 먼지 트랩 근처에서는 2.0 이상의 국지적 변동을 보인다.
  • NAUTILUS 모델과의 비교 결과, 관측된 높은 CS/SO 비를 재현하기 위해서는 C/O 비율이 $\gtrsim 1$인 모델만이 가능하다.
  • 데이터는 강한 황 결핍을 가리키며, 도출된 황 함량은 [S/H] $\approx 8 \times 10^{-8}$로, 우주 표준값($\sim 1.5 \times 10^{-5}$)보다 현저히 낮다.
  • 최적의 적합 모델 구성은 C/O 비율 1.0, 가스-먼지 질량비 40, 그리고 결핍된 황 함량을 포함한다. 그러나 저자들은 단일 모델이 관측된 모든 종(특히 H2CO 및 HCO+)의 방사형 프로파일을 동시에 성공적으로 재현하지 못했다는 점을 언급하였다.

의의 및 주장
본 논문은 AB Aur를 행성 형성 화학을 이해하기 위한 결정적인 사례 연구("로제타 스톤")로 위치시킨다. 주요 의의는 높은 CS/SO 비를 바탕으로 AB Aur 원반의 C/O 비율이 $\gtrsim 1$임을 도출해 냈다는 점에 있다. 이는 $\sim 93$ au에 위치한 후보 행성 AB Aur b를 포함하여 남은 원반에서 형성될 행성들이 태양계 행성들과는 구별되는 탄소가 풍부한 대기와 광물학을 가질 가능성이 높음을 의미한다.

저자들은 CS/SO 비가 C/O $\approx 1$과 하위 태양 수준의 황 함량을 지지하는 강력한 제약을 제공한다는 점을 인정하면서도, 현재의 화학 네트워크가 CS, SO, H2S의 함량을 동시에 재현하지 못하는 것은 황 함유 종의 형성 및 파괴 경로 중 중요한 부분이 현재 모델에서 누락되었음을 시사한다고 신중하게 결론지었다. 또한, 본 연구는 난류 및 비케플러 역학(예: 원반에 충돌하는 스트리머)이 분자 함량을 크게 변화시켜 정적인 원반 모델의 해석을 복잡하게 만들 수 있음을 강조한다.