핵붕괴 초신성 모델과 운석 동위 원소 데이터 비교를 위한 파이썬 툴 SIMPLE 소개

초록

본 논문은 핵붕괴 초신성(CCSN) 모델 18종을 포함한 데이터베이스와, 운석 및 CAI 등에서 측정된 동위 원소 변동성을 비교·시각화할 수 있는 파이썬 패키지 SIMPLE을 발표한다. 사용자는 모델별 동위 원소 질량분율을 손쉽게 불러와 몰분율로 변환하고, 태양계 표준과의 혼합선을 계산해 관측 데이터와 직접 대조할 수 있다. 예시로 Ni 동위 원소를 분석해, 운석 Ni 이례성이 CCSN 내부 Si‑rich 층에서 방출된 물질과 일치함을 확인한다.

상세 분석

SIMPLE은 현재 6개의 독립적인 연구팀이 발표한 15, 20, 25 M⊙ 질량의 핵붕괴 초신성 모델을 총 18개 수록하고 있다. 각 모델은 질량 좌표별(즉, 방출된 질량 구간) 동위 원소 질량분율(Xi)을 제공하며, SIMPLE은 이를 자동으로 몰분율(Yi = Xi/Ai)로 변환한다. 변환 과정에서 원자 질량을 정수(Ai)로 근사하고, 사용자 선택에 따라 질량분율↔몰분율 간 전환을 자유롭게 수행한다는 점이 핵심이다. 또한, 초기 별 조성(태양 금속성)과 모델별 초기 동위 원소 비율을 “표준”으로 삼아, 모델 자체 내에서의 상대적 변동만을 강조한다. 이는 실험실에서 사용되는 지구 표준과 혼동될 위험을 방지한다.

코드 구조는 크게 네 부분으로 나뉜다. 첫째, Zenodo에 저장된 텍스트/CSV 파일을 읽어 pandas DataFrame으로 변환한다. 둘째, 사용자가 지정한 원소·동위 원소 리스트에 대해 선택적 추출 및 결측치 경고를 제공한다. 셋째, 방출 질량 좌표를 x축으로, 선택된 동위 원소 비율(예: 60Ni/58Ni)을 y축으로 하는 혼합선 플롯을 자동 생성한다. 이때, “bulk Solar System” 조성(보통 지구 조성)과 모델 조성을 선형 혼합하여 다양한 희석 비율을 시뮬레이션한다. 넷째, 방사성 핵종 처리 옵션을 제공한다. 현재는 붕괴 시뮬레이션이 없으며, 방사성 핵종을 완전히 제외하거나 안정 핵종에 전부 합산하는 두 가지 선택만 가능하다.

SIMPLE의 장점은 (1) 모델별 동위 원소 프로파일을 직접 시각화함으로써 “층별 혼합”을 정량적으로 탐구할 수 있다. (2) 파이썬 기반이므로 사용자 정의 함수, MCMC 샘플링, 혹은 머신러닝 기반 파라미터 추정과 쉽게 연계 가능하다. (3) 오픈소스(GitHub)와 pip 배포를 통해 커뮤니티 기여를 장려한다는 점이다.

하지만 몇 가지 제한점도 존재한다. 현재 지원하는 모델은 모두 Fe‑core 붕괴 초신성이며, 전자 포획 초신성(ECSN)이나 3‑D 시뮬레이션 결과는 포함되지 않는다. 방사성 붕괴와 화학적 분별 효과를 정밀히 다루려면 별도의 후처리 모듈이 필요하다. 또한, 초기 조성 차이에 의해 발생하는 “인위적” 이례성을 완전히 제거하려면, 모든 모델을 동일한 초기 조성(예: Asplund 2009)으로 재계산해야 하는데, 이는 현재 데이터베이스에 포함되지 않았다.

핵심 과학적 인사이트는 Ni 동위 원소 분석 결과에서, 관측된 운석 Ni 이례성이 가장 안쪽 Si‑rich 층(보통 O‑Ne/Mg‑Si 혼합 구역)에서 방출된 물질과 일치한다는 점이다. 이는 Ni이 핵융합 단계에서 주로 Si‑burning 구역에서 생성되며, 해당 구역의 물질이 먼지 형태로 응축돼 원시 태양계 물질에 섞였을 가능성을 시사한다. 이러한 결론은 기존에 제시된 “단일 초신성 혼합” 모델을 정량적으로 뒷받침한다.

전반적으로 SIMPLE은 천체물리학자와 동위 원소 분석가가 모델과 데이터 사이의 격차를 메우는 실용적인 도구이며, 향후 모델 확대와 방사성 붕괴 모듈 추가를 통해 더욱 포괄적인 연구 플랫폼으로 성장할 잠재력을 가진다.