Constraining r-process nucleosynthesis with multi-objective Galactic chemical evolution models

기술적 요약: 다목적 은하 화학 진화 모델을 통한 r-과정 핵합성 제약

문제 제기
급속 중성자 포획 과정 (r-과정) 의 천체물리학적 기원은 여전히 불확실하며, 중성자별 병합 (NSMs) 과 자기-회전 초신성 (MR-SNe) 을 포함한 경쟁 시나리오들이 존재합니다. 다중신호 관측은 NSMs 를 유효한 발생 장소로 확인했지만, 그 지연된 발생과 낮은 빈도는 r-과정이 가장 초기의 저금속성 별에서 어떻게 풍부해졌는지를 설명하는 데 있어 NSMs 의 능력을 의문시합니다. 반면, MR-SNe 와 같은 즉각적인 원천은 이론적으로 타당하지만, 그 빈도와 수율에 대한 직접적인 관측적 제약이 부족합니다. 은하 화학 진화 (GCE) 모델링에서 중요한 쟁점은 'r-과정 보편성'이라는 가정, 즉 태양계 r-과정 분포를 기반으로 유로퓸 (Eu) 에서 다른 원소들까지 r-과정 원소의 상대적 풍부도 패턴이 균일하게 확장된다는 가설입니다. 저금속성 별에 대한 관측은 이 보편성이 무거운 원소 ($Z \ge 56$) 에 대해서는 성립할 수 있으나, 경량 중성자 포획 원소 ($38 < Z < 50$) 에 대해서는 실패할 수 있음을 시사합니다. 이러한 경량 원소들은 r-과정이 결핍된 별들에서 Eu 에 비해 종종 증가된 비율을 보입니다. 현재 모델들은 종종 단일 목적 최적화 (Eu 만 적합화) 에 의존하여, 여러 원소의 화학 진화를 동시에 재현하는 데 필요한 상충 관계를 가릴 수 있습니다.

방법론
저자들은 약 $1.5 \times 10^5$개의 단일 강하, 균질 은하 화학 진화 모델 그리드를 활용한 체계적이고 데이터 기반의 접근법을 사용합니다. r-과정을 특정 천체물리학적 장소 (예: NSMs 또는 MR-SNe) 로 고정하는 대신, 네 가지 주요 변수를 변화시키는 유연한 매개변수 프레임워크를 채택합니다:

1. 사건당 Eu 수율 ($Y_{Eu,r}$): $[10^{-7}, 10^{-5}] M_\odot$ 범위에서 로그적으로 샘플링.
2. 사건 비율 ($\alpha$): r-과정 물질을 생성하는 조상별의 비율로, $[0, 1]$ 범위에서 선형적으로 샘플링.
3. 지연 시간 ($\tau$): 항성 죽음과 풍부화 사이의 시간으로, $[0, 14]$ Gyr 범위에서 준로그적으로 샘플링.
4. 조성별 질량 범위 ($[M_l, M_u]$): $[8, 10)$에서 $[95, 100] M_\odot$까지 $5 M_\odot$ 간격으로 변화.

다른 중성자 포획 원소의 수율은 태양계 r-과정 패턴 (Simmerer et al. 2004) 에 따라 Eu 수율을 스케일링하여 도출됩니다. 모델들은 Gaia-ESO Survey, MINCE, SAGA 데이터베이스의 데이터를 결합하여 구성된 넓은 금속성 범위 ($-2.7 \le [\text{Fe/H}] \le \text{solar}$) 에 걸친 관측 데이터셋과 비교됩니다.

모델 성능을 평가하기 위해 저자들은 파레토 프론트 분석에 기반한 다목적 최적화 전략을 활용합니다. 단일 '최적 적합' 모델을 찾는 대신, 한 원소에 대한 적합도를 개선하면 반드시 다른 원소에 대한 적합도가 악화되는 파레토 최적 해의 집합을 식별합니다. 이를 통해 서로 다른 원소 제약 (예: Eu 대 Zr, 또는 경량 그룹 Sr-Y-Zr 대 중량 그룹 Ba-La-Ce) 간의 상충 관계를 매핑할 수 있습니다.

주요 기여 및 결과

• 매개변수 제약: 가장 성능이 우수한 모델들은 일관되게 짧은 지연 시간 ($\tau \le 30$ Myr), 저질량 조성별 (약 $20-25 M_\odot$ 중심), 그리고 사건당 $\sim 2 \times 10^{-7} M_\odot$의 유효 Eu 주입률을 선호합니다. $\sim 80 M_\odot$보다 무거운 조성별은 수율이 높더라도 관측된 r-과정 풍부화를 단독으로 설명하기에는 너무 드문 것으로 나타났습니다.
• 경량 원소에 대한 태양계 스케일링의 실패: 최적의 Eu 모델을 태양계 r-과정 패턴으로 스케일링할 때, 중간에서 저금속성 영역에서 경량 원소 (Sr, Zr) 의 풍부도를 체계적으로 과소평가하는 반면, 무거운 원소 (Ba, La, Ce) 는 더 성공적으로 재현합니다.
• 파레토 상충 관계:
  • Eu 대 Zr: 2 차원 최적화에서 단일 모델이 Eu 와 Zr 경향을 동시에 재현하지 못합니다. Zr 을 잘 적합하는 모델들은 체계적으로 Eu 를 과대평가하며, 그 반대도 마찬가지입니다.
  • 경량 그룹 대 중량 그룹: 경량 그룹 (Sr, Y, Zr) 에 대한 3 차원 최적화는 상당한 충돌을 드러내며, 파레토 프론트가 이상적인 해로부터 큰 거리를 보입니다. 반면, 중량 그룹 (Ba, La, Ce) 은 훨씬 낮은 상충 관계와 이상점까지의 더 작은 거리로 Eu 와 동시에 재현될 수 있습니다.
  • 전체 비교: Eu 를 경량 그룹 (Sr-Y-Zr) 과 결합할 때, 목적 함수 공간에서 이상점까지의 거리는 Eu 를 중량 그룹과 결합할 때 ($D \approx 50$) 보다 훨씬 크게 ($D \approx 200$) 나타나, 태양계 스케일링 관계가 경량 중성자 포획 원소로 갈수록 점차 유효성이 떨어진다는 것을 확인시켜 줍니다.

의의 및 결론
이 논문은 태양계 스케일 수율 가정 하에 작동하는 단일 클래스의 r-과정 사건으로는 은하 내 경량 및 중량 중성자 포획 원소의 화학 진화를 동시에 설명할 수 없다고 결론지었습니다. 관측된 긴장은 r-과정 보편성이 전체 질량 범위에 걸쳐 확장되지 않음을 시사합니다.

대신, 결과는 적어도 두 가지 다른 스케일링 구성 요소의 필요성을 시사합니다:

1. 주 r-과정 구성 요소 (H-구성 요소): 태양계 패턴과 유사하며 r-과정이 풍부한 별들에서 관측되는 것으로, 무거운 원소 ($Z \ge 56$) 의 생성을 지배합니다.
2. 약한 r-과정 구성 요소 (L-구성 요소): 란타노이드에 비해 경량 원소 (Sr, Y, Zr) 의 생성이 강화된 특징을 가지며, r-과정이 결핍된 별들에서 관측된 패턴과 유사합니다.

저자들은 이러한 구성 요소들이 서로 다른 천체물리학적 조건이나 사건에서 비롯될 가능성이 있다고 주장하지만, 구체적인 물리적 발생 장소를 특정하지는 않습니다. 이 연구는 화학 진화 모델에서 보편적인 r-과정 패턴을 가정하는 것의 한계를 강조하고, 복잡한 핵합성 기여를 분리해 내는 데 다목적 최적화의 유용성을 입증합니다. 향후 연구는 GCE 모델 내에서 약한 구성 요소를 직접 제약하거나 시뮬레이션된 풍부도 패턴을 통해 이를 규명할 것을 제안합니다.