태양 산소 풍부도와 태양 주기 변동 연구

초록

본 논문은 20년간의 Sun‑as‑a‑star 관측 데이터를 이용해 O I 적외선 삼중항(7772‑7775 Å)과 금지선(6300 Å)의 등가폭(EW) 변화를 태양 활동 지표인 태양반점수와 비교하였다. 고신호대비( S/N≈1000)인 PEPSI 데이터에서만 적외선 삼중항의 EW가 활동 최대와 최소 사이에 약 1 mÅ, 즉 산소 풍부도 A(O) 변동으로 환산하면 0.01 dex 정도의 미세한 증가를 보였으며, 이는 현재 문헌에서 보고되는 시스템적 불확실도보다 작다. 다른 장비(FEROS, HARPS‑N, NEID)에서는 통계적 유의미성을 찾지 못했다. 금지선 6300 Å에서는 어떠한 변동도 검출되지 않았다. 결론적으로, 현재 측정 정밀도 내에서는 태양 산소 풍부도가 태양 주기에 따라 실질적으로 변하지 않으며, 보다 활동적인 별에서는 이 효과가 편향 요인으로 작용할 가능성이 있다.

상세 분석

이 연구는 태양의 산소 풍부도 측정에 가장 널리 사용되는 O I 적외선 삼중항(7772, 7774, 7775 Å)과 비교적 약한 금지선(6300 Å)의 등가폭 변화를 장기적인 Sun‑as‑a‑star 스펙트럼 시계열을 통해 정량화했다. 사용된 네 개의 고해상도 분광기(FEROS, HARPS‑N, PEPSI, NEID)는 각각 2003‑2024년 사이의 데이터를 제공하며, 해상도는 R≈48 000–250 000, S/N은 200–1000 수준이다. 데이터 전처리 단계에서는 매월 15일 전후 3일 이내의 최고 S/N 관측을 선택하고, ARES v2 자동 EW 측정 코드를 이용해 가우시안 프로파일 피팅 및 로컬 연속선 정규화를 수행했다.

통계적 분석은 두 가지 접근법으로 진행되었다. 첫째, EW와 태양반점수(활동 지표) 사이의 피어슨 상관계수(r)와 p‑값을 계산하여 선형 상관성을 검증하였다. 둘째, 반점수 50을 기준으로 데이터를 ‘활동기(≥50)’와 ‘정지기(<50)’ 두 그룹으로 나누어 각각의 중앙값(EW_low, EW_high)을 구하고, 그 차이(ΔEW)를 산출하였다. ΔEW를 3D NLTE 모델(Steffen et al. 2015) 기반의 변환식에 적용해 산소 풍부도 차이(ΔA(O))로 환산하였다.

결과적으로, FEROS와 NEID에서는 r값이 0.1 이하이며 p‑값이 0.3 이상으로 귀무가설을 기각할 수 없었다. 반면 PEPSI 데이터는 r≈0.78–0.79, p < 0.001로 강한 양의 상관성을 보였으며, 활동기와 정지기 사이의 ΔEW가 0.8–1.0 mÅ, 이는 ΔA(O)≈0.01 dex에 해당한다. 이 변동은 현재 태양 산소 풍부도(A(O)≈8.69 dex) 측정에 사용되는 시스템적 오차(≈0.04–0.06 dex)보다 작다. 금지선 6300 Å는 PEPSI와 HARPS‑N 모두에서 r≈−0.06~0.16, ΔEW≈±0.03 mÅ 수준으로 통계적 유의미성을 찾지 못했으며, 이는 해당 라인이 활동에 민감하지 않음과 측정 정밀도 한계(≈0.05 mÅ) 때문으로 해석된다.

이러한 결과는 두 가지 중요한 함의를 가진다. 첫째, 현재 사용되는 고해상도 일일 관측 데이터와 전통적인 EW 측정 방법으로는 태양 자체의 산소 풍부도가 태양 주기에 따라 실질적으로 변한다는 결론을 내리기 어렵다. 둘째, 보다 활동적인 태양 유사 별(예: Hyades의 600 Myr 청색왜성)에서는 같은 메커니즘이 증폭되어 0.05 dex 이상, 혹은 그 이상의 풍부도 편향을 초래할 수 있다. 따라서 활동이 강한 별에 대한 산소 풍부도 추정 시, 적외선 삼중항 대신 금지선이나 다른 저활동 민감도 라인을 사용하거나, 활동 지표에 기반한 보정 모델을 적용하는 것이 필요하다.

추가적으로, PEPSI와 같은 초고신호대비 관측이 가능할 때만 미세한 EW 변동을 검출할 수 있음을 확인했다. 향후 연구에서는 NESSI와 같은 수치적 Sun‑as‑Star 통합 모델을 활용해 인위적인 활동 필링 팩터를 조절함으로써, 다양한 활동 수준에서의 스펙트럼 변화를 시뮬레이션하고, 관측적 보정 함수를 정밀하게 구축할 계획이다.