Charting circumstellar chemistry of carbon-rich asymptotic giant branch stars. III. SiO and SiS abundances

기술 요약: 탄소풍부 점성거성 (AGB) 항성의 성간 화학 지도 작성 III. SiO 및 SiS 풍부도

문제 및 배경
탄소풍부 (C 형) 점성거성 (AGB) 항성의 성간 껍질 (CSEs) 은 복잡한 분자 화학이 일어나는 장소입니다. 전형적인 탄소별인 IRC+10 216 은 광범위하게 연구되었으나, 탄소별 전체 집단 내 분자 풍부도에 대한 현재 이해는 주로 단일접시 (SD) 스펙트럼이나 제한된 간섭계 매핑에 의존하고 있습니다. 이러한 방법들은 종종 다양한 반경 범위에서 공간적 형태와 여기 조건을 동시에 제약하는 데 실패합니다. 특히, 규소를 포함한 난휘발성 종인 SiO 와 SiS 의 풍부도는 여전히 불확실합니다. C 풍부 환경에서 쉽게 형성되는 CS 와 달리, SiO 와 SiS 의 형성은 먼지 표면에서의 응축 (freeze-out) 이나 충격 화학과 같은 요인에 의존할 수 있습니다. 본 연구는 5 개의 탄소별 표본에 대해 SiO, SiS 및 주요 동위원소체 ($^{29}$SiO, $^{30}$SiO, $^{29}$SiS, $^{30}$SiS, Si$^{34}$S) 의 견고한 성간 풍부도 프로파일을 도출하여, IRC+10 216 이 전체 집단을 대표하는지 여부를 검증하고 화학 모델링 예측을 평가하는 것을 목표로 합니다.

방법론
저자들은 5 개의 고질량손실율 (MLR) C 형 AGB 항성 표본을 분석했습니다: IRAS 15194−5115, IRAS 15082−4808, IRAS 07454−7112, AFGL 3068, 그리고 IRC+10 216. 본 연구는 다음과 같은 포괄적인 데이터 세트를 결합하여 수행되었습니다:

• 공간적으로 분해된 데이터: ALMA 밴드 3 스펙트럼 탐사 (85–116 GHz) 와 IRC+10 216 의 경우 ALMA/IRAM 30m 온더플라이 (on-the-fly) 매핑을 결합한 데이터. 다양한 공간 규모에 걸쳐 선 플럭스를 추적하기 위해 컨볼루션된 큐브에서 다중 구경 스펙트럼을 추출했습니다.
• 단일접시 데이터: 낮은 J 에서 높은 J 전이까지 광범위한 여기 조건을 커버하기 위한 APEX, Herschel/HIFI, IRAM 30m, Yebes 40m 망원경의 스펙트럼.

저자들은 가속화 람다 반복 (ALI) 코드를 사용한 상세한 1 차원 국소 열평형 (non-LTE) 복사 전달 (RT) 모델링을 활용했습니다. 모델에는 다음이 포함되었습니다:

• 물리적 구조: 이전 먼지 및 CO 모델링 (Paper II) 에서 도출된 입력 매개변수 (항성 광도, 먼지 특성, 기체 밀도) 를 사용하여 일정한 MLR 과 등방성 팽창을 가진 구대칭 매끄러운 CSE 모델.
• 분자 데이터: SiO ($v=0, 1$ 포함) 와 SiS ($v=0, 1, 2$ 및 회전 - 진동 전이 포함) 에 대한 복사 및 충돌 여기율.
• 풍부도 프로파일: 피크 풍부도 ($f_0$) 와 e-감쇠 반경 ($R_e$) 을 자유 매개변수로 정의한 가우시안 프로파일.
• 제약 조건: 모델링된 선 적분 강도와 관측된 선 적분 강도 간의 $\chi^2$를 최소화하여 모델을 제약했으며, ALMA 다중 구경 스펙트럼의 공간적 제약과 SD 선의 여기 제약을 모두 활용했습니다.

주요 결과

1. 풍부도 프로파일: 공간적으로 분해된 데이터가 부족하여 $f_0$와 $R_e$값이 제약되지 않은 AFGL 3068 을 제외한 모든 소스에 대해 견고한 풍부도 프로파일이 도출되었습니다.

   • SiO: 4 개의 제약된 소스의 피크 풍부도는 $1.0 \times 10^{-7}$에서 $7.1 \times 10^{-7}$사이에 분포했습니다. 특히 IRAS 15194−5115, IRAS 15082−4808, IRAS 07454−7112 의 SiO 피크 풍부도는 IRC+10 216 ($2.2 \times 10^{-7}$) 보다 약 5 배 큰 것으로 나타났습니다.
   • SiS: 피크 풍부도는 표본 전체에서 더 균일하게 분포했습니다 ($2.0 \times 10^{-6}$에서 $4.7 \times 10^{-6}$). 이는 SiO 풍부도보다 일관되게 2~15 배 높았습니다.
   • 반경 범위: SiO 의 e-감쇠 반경 ($R_e$) 은 $1.3 \times 10^{16}$에서 $7.0 \times 10^{16}$cm 범위였으며, SiS 의 경우 $6.0 \times 10^{15}$에서 $1.0 \times 10^{17}$cm 범위였습니다. 두 종 모두 $R_e$가 기체 밀도가 증가함에 따라 증가하는 것으로 나타났습니다.

2. 동위원소 비율: 도출된 동위원소 비율 ($^{28}$Si/$^{29}$Si, $^{28}$Si/$^{30}$Si, $^{29}$Si/$^{30}$Si, $^{32}$S/$^{34}$S) 은 모든 소스에서 태양 값과 일치하여, 경미한 s-과정 증폭을 제외하고는 AGB 단계에서 이러한 원소들의 핵합성적 변화가 유의미하지 않음을 시사합니다.

3. 모델링의 어려움:

   • IRC+10 216 SiO: 저자들은 IRC+10 216 의 SiO 에 대해 낮은 J ALMA 선과 높은 J SD 선을 동시에 적합시킬 수 없었습니다. ALMA 데이터를 적합하는 모델은 SD 선을 과소평가했고, 그 반대의 경우도 마찬가지였습니다. 수정된 내부 풍부도 프로파일 (단계 함수) 이나 항성/먼지 매개변수 변동을 사용하여 이를 해결하려는 시도는 실패했으며, 이는 단순화된 항성 가정을 가진 1 차원 RT 모델링의 한계를 부각시켰습니다.
   • SiS J=6-5: SiO 와 유사하게, 모델은 SiS 의 $J=5-4$와 $J=6-5$선을 동시에 적합할 수 없었습니다. 저자들은 이를 $v=1$레벨에 영향을 미치는 다른 분자 (예: C$_2$H$_2$, HCN) 와의 IR 펌핑 중첩으로 인한 것으로 귀인하며, 이는 표준 RT 모델에 포함되지 않았습니다. 결과적으로 $J=6-5$선은 최종 풍부도 제약에서 제외되었습니다.

4. 화학 모델과의 비교:

   • SiO: 화학 모델은 도출된 SiO 풍부도 프로파일과 반경 범위를 잘 재현했습니다.
   • SiS: 화학 모델은 SiS 에 대해 $R_e$값을 일관되게 과대평가하여, 화학 모델에 사용된 현재 SiS 의 광분해율이 과소평가되었음을 시사합니다.

의의 및 결론
본 논문은 광범위한 여기 조건에 걸친 공간적으로 분해된 관측이 성간 분자 풍부도에 대한 견고한 제약을 필요로 함을 보여줍니다. 본 연구는 다른 소스들이 현저히 높은 SiO 풍부도를 보임에 따라 IRC+10 216 이 C 형 AGB 집단을 완전히 대표하지는 않는다는 것을 확인했습니다. 연구 결과는 SiO 와 SiS 선의 여기에서 적외선 펌핑의 결정적 역할을 강조합니다. 또한, 관측된 SiS 반경 범위와 화학 모델 예측 간의 불일치는 SiS 의 광분해율 수정이 필요함을 시사합니다. 저자들은 1 차원 RT 모델링이 효과적이지만 복잡하고 매끄럽지 않은 성간 환경에 적용될 때 내재된 한계가 있으며, 향후 AGB 화학에 대한 포괄적인 개요를 위해서는 질량손실율이 낮은 소스를 포함한 더 큰 표본과 업데이트된 화학 모델이 결합되어야 한다고 결론지었습니다.